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Fluctuation theorems (FTs), which describe some universal properties of nonequilibrium fluctu- 
ations, are examined from a quantum perspective and derived by introducing a two-point mea- 
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I. INTRODUCTION 

Small fluctuations of systems at equilibrium or weakly 
driven near equilibrium satisfy a universal relation 
known as the fluctuation-dissipation (FD) theorem 
[13, m, [mi, [HI [l67, 201]. This relation that connects 
spontaneous fluctuations to the linear response holds 
for classical and quantum systems alike. The search for 
similar relations for systems driven far from equilibrium 
has been an active area of research for many decades. 
A major breakthrough in this regard had taken place 
over the past fifteen years with the discovery of exact 
fluctuation relations which hold for classical systems far 
from equilibrium. These are collectively referred to as 
fluctuation theorems (FTs). In order to introduce these 
theorems we will adopt the following terminology. A 
system that follows a Hamiltonian dynamics is called 
isolated. By default, we assume that the Hamiltonian is 
time independent. Otherwise, it means that some work 
is performed on the system and we denote it driven 
isolated system. A system that can only exchange 
energy with a reservoir will be denoted closed. If parti- 
cles are exchanged as well, we say that the system is open. 

The first class of FTs (and the earliest discovered) 
deal with irreversible work fluctuations in isolated 
driven systems described by an Hamiltonian dy- 
namics where the Hamiltonian is time- dependen t 

[11 [Ti, fil liO, [M [si, [S [13, [M [ii, [al [9^^ [Tm ITesl. 



An example is the Crooks relation which states that the 
noncquilibrium probability p(W), that a certain work 
w = W is performed by an external time-dependent 
driving force acting on a system initially at equilibrium 
with temperature divided by the probability 

p{—W), that a work w — —W is performed by the 
time-reversed external driving force acting on the 
system which is again initially at equilibrium, satisfies 
piW)/p{-W) = exp[P{W - AF)], where AF is the 
free energy difference between the initial (no driving 
force) and final (finite driving force) equilibrium state. 
The Jarzynski relation (exp [— = exp [—(3AF] 
follows immediately from J dWp{—W) = 1. A second 
class of FTs is concerned with entropy fluctuations in 
closed systems described by deterministic thermost atted 
equations of motions [H [H [H, [H, [zl, [tI, \iEt\ 
and a third class treats the fluctuations of entropy (or 
related quantities such as irreversible work, heat and 
matter currents) in closed or open systems described 
by a stochastic dynamics [IJsijll EH, [H, El, [H, 
[73, M, [HE [m [iM fl60l.TlMll77f. As an example 
for the last two classes, we give the steady-state FT 
for the entropy production. We consider a trajectory 
quantity s whose ensemble average (s) can be associated 
with an entropy production (the specific form of s 
depends on the underlying dynamics). If p{S) denotes 
the probability that s = S when the system is in a 
noncquilibrium steady-state, then for long times the FT 
reads p{S)/p{—S) = exp [S]. FTs valid at any time such 



as the work FTs are called transient FTs while those 
who require a long time limit are called steady-state FTs. 

The FTs are all intimately connected to time-reversal 
symmetry and the relations between probabilities of 
forward and backward classical trajectories. Close to 
equilibrium the FTs reduce to the known fiuctuation- 
dissipation relations such as the Green- Kubo relation 
for transport coefficients [1, i, [fii, [fii, [ul- These 
classical fluctuation rela tions have been reviewed in 
Refs. tZCi JX, IE, M [Hi]. Some of these relations were 
verified experimentally in mesoscopic systems where 
fluctuations are sufficiently large to be measurable. 
Work fluctuations have been studied in macromolecule 
pulling experi ment s (40l . Il24l | and in optically driven 
microspheres jl82j |. entropy fluctuations have also 
been measured in a similar system [l88j and in spec- 
troscopic experiments on a defect center in diamond 
[I59I Il79j |. When decreasing system sizes, quantum 
effects may become significant. Applying the standard 
trajectory-based derivations of FTs to quantum regime 
is complicated by the lack of a classical trajectory 
picture when coherences are taken into account and by 
the essential role of measurements, which can be safely 
ignored in ideal classical systems. We show that the FTs 
follow from fundamental dynamical symmetries that 
apply equally to classical and quantum systems. 

Earlier derivations of the Jarzynski relation were 
done for quantum systems b y de f ining a w ork operator 
[E M, MM, MMMM till [TM lilg . since work 
is not in general an ordinary quantum "observable" (the 
final Hamiltonian does not commute with the initial 
Hamiltonian) jl76l |. attempts to define such an operator 
had led to quantum corrections to the classical Jarzynski 
result. However, the Jarzynski relation in a closed driven 
quantum system may be derived without quantum cor- 
rections by introducing an initial and final projective 
measurment of the system energy in accordance with 
the quantum mechanical measurement postulate. This 
has been done (no t always in a e xplicit way) in Refs. 
[nel fml . HM flTi . ITtI Irzel . llTSt . The work is then a 
two-point quantity obtained by calculating the difference 
between the initial and final energy of the system. 
When the reservoir is explicitly taken into account, 
the Jarzynski relation has often bee n de rived using a 
master equation approach [44, 4^, 56. llSOll . Alternative 
derivations can be found in Refs. n3U]l73l|. 



The derivation of a steady-state FT for 
quantum s yste ms has been considered as well 
[i,i,[3SESM[lS[8i,[9i,[l4i,[M[lll[Ill[l83. Be- 
cause of the need to describe nonequilibrium ffuctuations 
in closed or open quantum systems exchanging energy or 
matter with their reservoir, many similarities exist with 
the rapidly dev eloping field of electroi i counting stat i stics 
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Il83l . Il86l Il9l| , where small nano-scale electronic devices 
exchange electrons. Fluctuations in such systems can 
nowadays be ex p erimen tally resolved at the single elec- 
tron level {2^, IGo, [13, [13, 1125 *1 . Similarities also exist with 
the more established field of photon counting statistics 
where photons emitted by a molecule or an atom driven 
out of equilibrium by a laser, arc individually detected 

[n,, 105,, IM,'a27„ 128,, OM, 1I55, JM, .198.] . 

Different types of approaches have been used to derive 
these FTs and describe these counting experiments. The 
first is based on the quantum master e quation (QME ) 

[m, m, y, m, m, m, m, m [m [m ee^ [m, tip. 

Here one starts with an isolated system containing 
the system and the reservoir in weak interaction. By 
tracing the reservoir degrees of freedom, taking the 
infinite reservoir limit and using perturbation theory, 
one can derive a closed evolution equation for the 
reduced density matrix of the system. The information 
about the reservoir evolution is discarded. However, 
the evolution of a quantum system described by a 
QME can be seen as resulting from a continuous 
projective measurement on the reservoir leading to a 
continuous positive operator-valued measurement on 
the system. Such interpretation allows to construct a 
trajectory picture of the system dynamics, where each 
realization of the continuous measurement leads to a 
given system trajectory (23, 24, 25, 76, 141]. The QME 
is recovered by ensemble averaging over all possible 
trajectories. This unraveling of the QME into trajec- 
tories has been originally developed in the description 
of photon counting statistics [tI [m^, [IM fl9i |. 
Another approach is based on a modified propaga- 
tor defined on a Keldysh loop which, under certain 
circumstances, can be interpreted as the generating 
function of t he electron coun t ing probab ility distribution 
[H, [H, M, [103, \lM, [133, [iM [Hi. Using a path 
integral formalism, the propagator of the density matrix 
of a "detector" with Hamiltonian /2m interacting 
with a system, can be expressed in term of the influence 
functional that only depends on the system degrees of 
freedom [b^]. The modified propagator is the influence 
functional when the system is linearly coupled to the 
detector (with coupling term xA, where x is the position 
of the detector and A a system observable) in the limit 
of very large detector inertia to — > 00. It is only under 
some specific assumptions (such as a classical detector 
where the detector density matrix is assumed diagonal) 
that the modified propagator becomes the generating 
function associated with the probability distribution that 
the detector momentum changes from a given amount 
which can be interpreted as the probability to measure 
the time average of the system observable A: drAlr). 
If A is an electric current, then the integral gives the 
number of electrons transfered. An early quantum FT 
for electro nic j unctions has been derived in this context 
in Ref. [l80| based on the time-reversal invariance 
of the Hamiltonian quantum dynamics. Different 



derivations of quantum FTs relying o n this approach 
have been considered in Ref. jl53l Il54| . A third, semi- 
classical scattering, approach is often used in electron 
counting statistics [ToO, [HE [HI [Til [TUl. This can 
be recovered from the modified propagator approach 
as recently shown in |164| |. but will not be addressed here. 

We consider fluctuations in the outputof a two-point 
projective measurement (of energy, particle, charge, etc.). 
This allows us to avoid the detailed modeling of detectors 
and their dynamics. The projective measurement can 
be viewed as an effective modeling of the effect of the 
system-detector interaction on the system or as resulting 
in a fundamental way from the quantum measurement 
postulate. The three other approaches (unraveling of 
the QME, modifled propagator on Keldysh loop and the 
scattering approach) can be recovered in some limits of 
the two-point measurement approach. This provides a 
unifled framework from which the different types of FTs 
previously derived for quantum systems can be obtained. 

In section [III we give the general expression for the 
probability of the output of a two-point measurement at 
different times on a quantum system described by the 
quantum Liouville equation. The calculation is repeated 
for a system described by the time-reversed dynamics. 
In section IIIIl we start by discussing the basic ingredi- 
ents required for FTs to hold. We use these results to 
derive three transient FTs, the Jarzynski and Crooks re- 
lation in isolated and closed driven systems and a FT 
for matter and heat exchange between two systems in di- 
rect contact. We also show that a steady-state FT can 
be derived for matter and heat exchange between two 
reservoirs through an embedded system. In section IIVI 
we consider a small quantum system weakly interacting 
with multiple reservoirs. We develop a projection super- 
operator formalism to derive equations of motion for the 
generating function associated with the system reduced 
density matrix conditional of the output of a two-point 
measurement of the energy or number of particles in the 
reservoirs. We apply this generalized quantum master 
equation (GQME) formalism to calculate the statistics 
of particles or heat transfer in different models of general 
interest in nanosciences in order to verify the validity of 
the steady-state FT. In section |Vl we present a nonequi- 
librium Green's functions formalism in Liouville space 
which provides a powerful tool to calculate the particle 
statistics of many body quantum systems. In section 
IVIl we show that the FTs can be used to derive gener- 
alized fluctuation-dissipation relations. Conclusions and 
perspectives will be drawn in section fVIII 



II. TWO-POINT MEASUREMENT STATISTICS 

We consider an isolated, possibly driven, quantum sys- 
tem described by a density matrix p{t) which obeys the 
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von Neumann (quantum Liouville) equation 

Its formal solution reads 

p{t)^U{t,Q)poU\t,{)) . 
The propagator 



U{t,Q) = exp+{-- dTH{T)} 

n=l "^0 "^0 



(1) 



(2) 



(3) 



dtr, 



H{h)H{t2)...H{tn) . 

is unitary {t,Q) = U-^{t,Q) and satisfies U\t,Q) = 
U{0,t) and U{t,ti)U{ti,0) = U{t,0). We use the sub- 
script + (— ) to denotes a antichronological (chronolog- 
ical) time ordering from left to right. We call ^ the 
forward evolution to distinguish it from the the time- 
reversed evolution that will be defined below. 



A. The forward probability 

We consider an observable A{t) in the Schrodinger 
picture whose explicit time dependence solely comes 
from an external driving. For non-driven systems 
A(t) = A. In the applications considered below, A{t) 
will be either an energy operator H or & particle number 
operator N . The eigenvalues (eigenvectors) of A[t) are 
denoted by at (|at)): A{t) = Y^a, Wt)at{at\- 

The basic quantity in the following discussion will be 
the joint probability to measure ao at time and at at 
time t 

P[at,a^] = Tr 

= P*[at,ao] , (4) 



where the projection operators are given by 
Pat = \at){at \ ■ 



(5) 



Using the properties Pa, = P^, and J2at ^o-t = 1, we 
can verify the normalization X^atao ^t'**' °'o] ~ C!on- 
sider two complete Hilbert space basis sets {|i,ao)} and 
{|jiflt)}i where i (j) are used to differentiate between 
the states with same Oq (at). The basis {|i,ao)} is cho- 
sen such that it diagonalizes /5o (this is always possible 
since po is hermitian). We can also write (jH) as 



at] I, floj 



(6) 



where 

P[j,at;i,aa] = | (j, at |[/(t, 0)|i, ao)P(i, ao|po|«, ao) • (7) 

The probability distribution for the difference Aa — 
at — ao between the output of the two measurements is 
given by 



p(Aa) = ^ (5(Aa - (a* - ao))P[af , ao] , 



(8) 



where 6{a) denotes the Dirac distribution. It is of- 
ten more convenient to calculate the generating function 
(GF) associated with this probability 



G(A) 



dAa 



.iAAa 



p{Aa) = G*(-A) 



J2e'^^-'-^"^P[at,ao] 



(9) 



The n'th moment, (Aa"), of p{Aa) is obtained by taking 
7i'th derivative of the GF with respect to A evaluated at 
A = 



(Aa") 



gn 

(_i)« -^GiX) 



A=0 



We further define the cumulant GF 
Z(A) = InG(A) . 



(10) 



(11) 



The n'th cumulant, of p(Aa) is obtained by taking 
n'th derivative of the cumulant GF with respect to A 
evaluated at A = 



Kn = (-i)" ^Z(A) 



(12) 



A=0 



The first cumulant coincides with the first moment which 
gives the average Ki — (Aa). Higher order cumulants 
can be expressed in term of the moments. The variance, 
K2 — (Aa^) — (Aa)^, gives the fluctuations around the 
average, and the skewness if 3 — ((Aa — (Aa))'^ gives 
the leading order deviation of p(Aa) from a Gaussian. 
When measuring the statistics of quantities associated to 
nonequilibrium fluxes, in most cases (but not always [55]) 
the cumulants grow linearly with time and it becomes 
convenient to define the long time limit of the cumulant 
GF 



5(A) = lim -Z{\) 

t~*oc t 



(13) 



whic h me asures the deviations to the central limit theo- 
rem [165| . 

We next turn to computing the GF. The initial density 
matrix can be expressed as 



Pa = Pq + Pa 



where 



Po 



1,3 



2^ PaoPoPa'^ 



(14) 



(15) 
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po commutes with A{0). Using the fact that f{A) = 
Pafio) where / is an arbitrary function, and using 

also 

Ee"''"°^<^oPoPao = e-'*-^(°)p^oe-'*^(°) , (16) 

we find, by substituting Q in lU, that 

G(A) = Trp(A,0, (17) 
where we have defined 

p{\t) = U^{t,0)~paUl,{t,Q) (18) 

2 2 

and the modified evolution operator 

U),{t, 0) = e'^^(*)c/(t, 0)e-'^^(") . (19) 

For A = 0, p{X,t) reduces to the system density matrix 
and U\{t, 0) to the standard evolution operator. Defining 
the modified Hamiltonian 

Hx{t) = e'^^(*)i?(Oe-'^^(*^ - hXdtA{t) , (20) 

we find that U\{t, 0) satisfies the equation of motion 

4f>A(<,0) = -^i?A(t)f/A(i,0). (21) 
at n 

Since U\{0,0) = 1, we get 

Ux(t,0) = exp, {-^ / dTHx(T)} (22) 
" n Jo ^ 

i /■* 

uK{t,0) = exp_{- / dTH_x{T)} . (23) 

2 n Jo ^ 

Equations ^ and ^ together with ^ and ^ 
provide an exact formal expression for the statistics of 
changes in A{t) derived from the two-point measure- 
ments. 

We note that if and only if the eigenvalues of A are 
integers (as in electron counting where one considers the 
number operator), using the integral representation of 
the Kronecker Delta 

can be written as 

P(A,0= / ^ p(A,A,t), (25) 

where 

p(A, A,t) = U^^x{t,Q)poU{_,{t,Q) . (26) 

2 ^^2 

We see that by introducing an additional A dependence, 
we where able to keep the initial density matrix po in 
(I^S)) instead of po as in (fT5|) . 



The current operator associated with A{t) is given by 
i{t)^'-[H{t),A{t)]+dtA{t). (27) 

As a result, 

i^'^Xt) = j^A^^\t) , (28) 

where the subscript [h) denotes the Heisenberg represen- 
tation i('*)(i) = ij^{t,{))A{t)lJ{t,Q). We can write ^ 
as 

Hx{t) = H{t) - \ni{t) + 0{\^h^) . (29) 

In the semiclassical approximation where terms 0{\^fi?) 
are disregarded, the GF ^ [with HI]), ^ and 
after going to the interaction representation, becomes 

G(A) = Tr{eJ ^0 '^^^<'''M^v;#/o -^^^''-'M} . (30) 

This form is commonly found in the modified propagator 
approach (described in the introduction) to counting 
statistics [iOTI. [lOSllTsgt . Notice that in these Refs. the 
full initial density matrix po is used in (j30p instead of 

Po- 

In most applications considered in this review, we will 
consider initial density matrices with no initial coher- 
ences in A{Q) space 

[i(0),po]=0. (31) 

This is equivalent to say that [Pao , Po] =0 or that po = 
Pq. In this case, Eq. Q can be written as 

P[at,ao]=Tr\^U\t,0)Pa,U{t,0)Pa„po} (32) 
and using ([5^ in 0, the GF simplifies to 

G(A) = Tr J^e'^uHt,o)Ait)uit,o)^-ixAio)~^^ (33) 

B. The time-reversed probability 

The time-reversed evolution brings the final density 
matrix of the forward quantum evolution ^ back to 
its initial density matrix. This means that if the initial 
condition of the time-reversed evolution is p^ — p{t) — 
U{t, Q)poU^{t, 0), the time-reversed evolution must be de- 
fined as p^'it) = C/t(i,0)p*o'C/(t,0), so that p^'{t) = po- 
The time-reversed expression of the two-point probability 
(|4]) is therefore 

P''[ao,at]=TT{Pa„U\t,0)PaJl'Pa,Uit,0)Pa,} .(34) 

A more systematic discussion on time-reversal operation 
in quantum mechanics and its relation to the definition 
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(|34)) is given in appendix]^ Without loss of generality, 
we choose a basis set {\j,at)} that diagonalizes /5q'', to 
show that (l34l) can be rewritten as 



ao;j,atJ , 



(35) 



where 

P*'-[i,ao;j,at] = \{j,at\U{t,0)\i,ao)\^{j,at\p'J\j,at) .(36) 

The probability to measure the difference Aa = a^ — at 
between the two measurements is given by 

p'^(Aa) = ^<5(Aa-(ao-at))P*^[ao,at] . (37) 
The associated GF reads 



(38) 



G^'{\) = / dAa e'-^^y'^(Aa) 

J — OQ 

= ^ e-'-^('^*-'^«'P"-[ao,at] 



Let us note that for a non-driven system with — po, 
using (|4]) and (p4)) . we find that P[at,ao] = P*''[at,ao]. 
This means, using ^ and ([35]). that 

(39) 
(40) 
(41) 



and 



p*'(Aa) =p(Aa) 



G*^(A) = G(A) 



Using again the partitioning 



Po=~Po+~Po 



where 



at at^a\ 

and following the same procedure as for the forward GF, 
we obtain 



G'^(A) =Trp"(A,t) 



where 



p*'-(A,t) = C/l(i,0)^^'-C/_^(t,0) 



(43) 
(44) 



As for in the semiclassical limit wc find 

G-(A) = Tr{eL^^°'^^''''^^^^!i^e/^" '^^^""(^)} . (45) 

We again note that if the initial density matrix of the 
time-reversed evolution contains no initial coherences in 
A{€) space 

[i(i),P^'-]=0, (46) 
or equivalently if [Pat,Po ] = or p^^ = Pq'', (I34p becomes 

P*'[ao,at] = Tr{c/(t,0)P„„?7^(t,0)Pa,Po } , (47) 



and 



G 



*'XA) = Tr |ei^^>(*'0)^(o)^>t(*.o)e-'^^(*)p*'j . (48) 



III. THE FLUCTUATION THEOREM 
A. General derivation and connection to entropy 

We define the log of the ratio of the forward and time- 
reversed probabilities defined in section |TT1 which in the 
classical theory of FTs is associated with the irreversible 
contribution to an entropy change 



It follows from ^ and ^ that 

(i, ao|/5o|i, ao) 



0', at|Po'lj>t> 



(49) 



(50) 



An integral FT immediately follows from the normaliza- 
tion of P*'' [i, ao; j, at] 

j,at.i,ao 

= ^*'[i,ao;j- at] = l. (51) 

j,at,i,ao 

Using Jensen's inequality (e"^) > e^^\ ([?T|) implies 

(P) = Y ^t^' "o] ^t^' ao] > . (52) 

j,at,i,ao 

Using we see th at (R) res embles a Kullback-Leibler 
(or relative) entropy [ll3l Il4l| . 

We define the probability distributions 

P{R)= Y Pij'^t''i''^o]HR- R[j,aui,ao]) (53) 

j,at,i,ao 

p*'(P)= J2 P''[hao;j,cit]6{R~R''[t,ao;j,at]){M) 



where 



P'''[i,ao; j, at] = In 



(55) 



(56) 



P[j,at;i,ao\ 
Using ([50]) and ([55)) . we see that 

P"'[z,ao; j, at] = -P[j, at; i, ao] . 
It then follows that 
p{R) = 

Y e^[^'"'^^''^''lp''^[z, ao; J, at]S{R - P[j, a*; z, ao]) 



],at,i,ao 



gi? ^ P*''[i,ao; j,at](5(P- P[j,at;?,ao]) 

J, at ,i,oo 

gfl ^ pt^[z,ao;i,at]<5(P + P*^[i,ao;j, at]) 

j,at,i,aQ 

eV(-P) , (57) 
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which gives the detailed FT 
p{R) 



In 



pt^(-i?) 



= R 



(58) 



The FTs (I5l1) and ((58)) are completely general but 
only useful when R can be exclusively expressed in 
terms of physical and measurable quantities (the 
eigenvalues of A{0) and A(t)). In sections IIII Bl and 
IIII C[ we will see that the i and j dependence of R, 
that labels states which cannot be differentiated by a 
projective measurement of the physical observable A{t), 
can be eliminated by making specific choices of p and p'"'. 

If the assumptions ([3T|l and (j46| are satisfied (this will 
be the case in most of the following applications), ([5^ 
can be expressed in term of quantum entropies. Using 
([50)) . the general property 



at; I, floj 



{i,ao\po\i,aa) 



(59) 



and the fact that [using assumption (PT|) ] 

^P[j>t;i,ao] = {j,at\p{t)\j,at) , (60) 



(I52p can be rewritten as a quantum relative entropy j23l . 
Il4l| between p{t) and p^J 



{R) =S-S = Ti-p{t){\np{t)-\npl') >0, (61) 



where 



and 



-Trp(t) ln/5(t) = -Trpolnpo 

- ^(i, ao\po\i, ao) ln(i, ao|/5o|z, ao) (62) 



-Tr/5(t) ln/5^'-==-Tr/5oln/5*^(t) 



,at\pl'\j,at) . (63) 



The second line of (|62p [ ([55]) ] is obtained using the as- 
sumption pip [ l|46p ]. S" is a von Neumann entropy but S 
is not. It can be compared to the von Neumann entropy 

S'^- EE -Trp''-(t)ln/5*"(t) = -Trp*^lnp*^ 

= -J2^J,at\Po\j,at)\n{j,at\pl'\j,at) (64) 

which is obtained using the general property 

P"'[i, ao;j, at] = {j, atlPo'li, at) (65) 

i,ao 

together with [using assumption 

YP*-''[i,ao;j,at] = {i,ao\p'' {t)\i,ao) . (66) 



We will see in the following applications that (i?) is 
always associated to the irreversible contribution of an 
entropy change. Eq. ([6T|l is therefore the qua i itum 
analog of the classical relation derived in Refs. [s^ . Il04l | 
and of the stochastic relation of Refs. [tI Fzst. 

In appendix m following Refs. [28l.[l95l|. we show that 
if one allows for a coarse-graining of po and pg in their 
measured subspaces, one can derive FTs for i?'s which 
can be expressed exclusively in terms of measurable 
probabilities (no i and j index) and such that (R) is 
the difference between the Gibbs-von Neumann entropy 
associated to the coarse-grained pg'' and /5g. 

We now examine the detailed FT from the GF per- 
spective. We define the GFs associated with p{R) and 
p*'(i?) 

/OO 
dR e'^%{R) 
-OO 

/OO 
dR e'^V''(^) • 
-OO 



(67) 



By combining ([55]) with (p7)) . we get 

G(A) = G"^(i- A) . 



(68) 



For a non-driven system with p^J = pg, we have seen that 
(1301) is satisfied. Combining this with the detailed 
FT ((58)) implies the fundamental symmetry G(A) — G(i— 
A) on the GF. This type of symmetry will be used in 
section IVII to derive generalized fluctuation-dissipation 
relations. 



B. Transient fluctuation theorems 



In this section, we show that the FT ([58)) can be used 
to derive the Crooks [i^, |4^ [o^l and the Jarzynski rela- 
tions [9^, in either isolated or closed driven quan- 
tum systems as well as a FT for for heat and particles 
exchange between two flnite systems. 



1. Work fluctuation theorem for isolated driven systems 

We consider an isolated system initially described by 
the Hamiltonian H{0) and at equilibrium e~''-^*^°^/Z(0), 
where Z{0) = Tre"'^'^'-"'' is the partition function. We 
can imagine that the system was in contact with a reser- 
voir at temperature /3~^ for t < 0. At t — the reser- 
voir is removed and the system energy is measured for 
the first time. After the first measurement, the system 
is then subjected to an external and arbitrary driving 
(the Hamiltonian is time-dependent). The second energy 
measurement occurs at time i, where the Hamiltonian is 
H(t). From the two measurements of this forward pro- 
cess we can calculate P[Et, Eq]. 
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In the backward process, the isolated system is initially 
described by the Hamiltonian H{t) and at equilibrium 
g-/5if(t)/2-(t), where Z{t) = Tre"'^^^*). We can imagine 
that at the end of the forward process, the system de- 
scribed by the Hamiltonian H{t) is put in contact with 
a reservoir at temperature f3~^ until it thermalizes, and 
that the reservoir is then removed at time zero when the 
energy of the system is measured for the first time in the 
backward process. After this first measurement, an ex- 
ternal driving, which is the time reversed driving of the 
forward process, is applied. The second energy measure- 
ment occurs at time t, where the Hamiltonian is H{0). 

In appendix |3 we show that the time-reversed evo- 
lution (as defined in section III Bp of an isolated system 
driven externally according to a given protocol, corre- 
sponds to the forward evolution of the isolated system 
externally driven according to the time-reversed proto- 
col. This means that the backward process just described 
is identical to the time-reversal of our forward process, 
so that the two measurements occurring during the back- 
ward process can be used to calculate P'''[£'o, Et]. 

To make the connection with the results of section llll 
we define the initial density matrices for the forward and 
backward process 



Po 



PO 



(69) 



We further set at — Et (ao — Eq), where H(t)\Et,j) — 
Et\Et,j) {H{0)\Eo,i) = Eo\Eo,i)). The index j (i) dis- 
tinguish between degenerate eigenstates so that {\Et,j)} 
{{\EQ,i)}) constitute a complete basis in Hilbert space. 
We also define the free-energy difference AF(t) = F{t) — 
F(0) between the initial and final state, where F(t) = 
—l3~^\nZ{t). Since the system is isolated, no heat ex- 
change occurs and the change in the system energy can 
be interpreted as the work done by the driving force on 
the system 



w = Aa = Et — Eq 
Eq. and ^ become 



(70) 



P[j,Et;i,Eo] = \{j,Et\U{t,0)\hEo)\'eM^"-^^°^) 

P'^[i,Eo;j,Et] = \{j,Et\Uit,0)\i,Eo)\'e-^(''^-^^'^) , 

so that Eq. ([50)1 becomes 

R[j, Et; I, E^] = (3{w ~ AF{t)) = R[Et, Eq] . (71) 

The essential property that R is independent of i and j 
and only expressed in terms of observable quantities is 
therefore satisfied. 

and become 



S = P{TrHit)pit) - F{t)) 
S = f3{TrH{0)po - F(0)) , 



and 



where 



{R)=S-S = /3((w) AF) > 



TvHit)p{t) - TTHiO)pa 



(73) 



(74) 



(w) is the average work, so that /3^^(i?) is the irre- 
versible work (the irreversible contribution to the entropy 
change). Using we get the Crooks relation 



Pjw) 
p*''(— w) 



AF) 



(75) 



The Jarzynski relation follows immediately from (j75p [by 
integrating p*''(— w) over w which is equal to one because 
of normalization! 



-!3w\ 



-I3AF 



(76) 



Eq uations ([75)1 and ([75]) have been first derived in Ref. 

1161 for a periodic driving (where AF = 0) and in Ref. 

178l | for finit e AF. Furth er studies of (l76|) ha ve b een 
done in Refs. [Ulllzlllzi and of ^ in Refs. ^7^. It 
was gen eralized to the microcanonical ensemble in Refs. 
f29l.[l75l. 



2. Work fluctuation theorem for closed driven systems 

We consider the same forward and backward process 
as described above, except that during the driving the 
system now remains in weak contact with a reservoir at 
equilibrium. The total Hamiltonian is therefore of the 
form H{t) = Hsit) + Hb + V, where Hs{t) {Hb) is the 
system (reservoir) Hamiltonian and V the weak interac- 
tion between the two. The work done by the driving 
force on the system is now given by the difference be- 
tween the system and the reservoir energy change (this 
last one represents heat) according to the first law of 
thermodynamics. 

In this case, the connection with the results of section 
in] is done using 



Po 



ZsiO) 



Po = \ \ ■ (77) 



Zs{t) Zb 



as well as ao — Eg (0) -I- Et and at — Es> (t) + Ef , where 
Es{0) {Es{t)) are the eigenvalues of HsiO) {Hs{t)) and 
Eb the eigenvalues of Hb- We define i = {is, it) and j = 
{js,jb), where is and js are used to distinguish between 
degenerate eigenstates of J?s(0) and Hs{t) and and 
jb between degenerate eigenstates of Hb- The work is 
therefore 



Aa = Us's + Q. 



b'b 



(78) 



(72) 



where Us's = Es'{t) — Es{0) is the change in the system 
energy and Qb'b — Eb> — Eb is the heat transferred from 
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the system to the reservoir. Since the eigenstates of the 
Hamiltonian constitute a complete basis set, ([7]) and ([55]) 
become 

P[j,E,,{t)+Ei,r,i,EM+Eb] (79) 
^\{js'b'\U{t,0)\tsb)\^{sb\po\sb) 

P''[i,EM + Eb;j,E,,it)+Eb^] (80) 
^\{Js'b'\U{t,0)\^sb)f{s'b'\pl^■\s'b'). 

Eq. ([50|l therefore gives 

R[j, Es. [t) + Eb' ; i, Es{0) + Eb] = f3{w - AF) 

= R[E,,{t)+Eb',E,{0) + Eb],m 

where AF(t) — F{t) ~ F{0) is the free-energy differ- 
ence between the initial and final system state {F{t) — 
— /3~^ In Zs{t)). The essential property that R is indepen- 
dent of i and j and expressed solely in terms of observ- 
able quantities is therefore again satisfied. Using ([58]) . 
we get the same Crooks (|75p and Jarzynski (|76p relation 
as in the isolated case. The two relations were derived 
for quantum open d riven systems in ma ny different ways 
in Refs. [3 EE S [Til EHl [US, Eli ■ using §^ and 
we also find that (USD still holds with 

{w) = TT{Hs{t) + HB)m - Tr(7?s(0) + Hb)po ■ (82) 



3. Fluctuation theorem for direct heat and matter exchange 
between two systems 

We consider two finite systems A and B with Hamilto- 
nians Ha and Hb, each initially at equilibrium with its 
own temperature and chemical potential. The two sys- 
tems are weakly interacting, allowing heat and matter 
exchange between them. The total Hamiltonian is of the 
form Htot — Ha + Hb + V, where V is the coupling term 
between A and B. The joint Hilbert space is Ha x "Hb- 
The energy Ea and the number of particles ua of sys- 
tem A is measured at time zero and again at time t. We 
assume 



Po = Po = p2{Pa,Pa)Pb{Pb,Pb) 



where 



-px(Hx-lJ.xNx) j-^ 



X 



(83) 



(84) 



and X = A, B. Sx is the grand canonical partition 
function. The index ix is used to distinguish between 
eigenstates of Hx with same energy Ex and number 
of particles nx- We define i = (iaj^b) and a = 
{EA,nA,EB,nB)- Using (O and ([36|), we find 

P[i,a]i,a\ = |(i',Q;'|C/t|z,a)|^(Q;|/5o|a) (85) 
P"-[z,a;z',a'] = ,a'\Ut\i,a)\^{a'\poW) . (86) 

Eq. (inni) with 

R[a',a] = 



-Pa^Ea 
-Pb{{Eb 



E'a) ~ fiAiuA - n'A)) (87) 



Conservation laws imply that changes in matter and en- 
ergy in one system are accompanied by the opposite 
changes in the other system so that 



Ea- 

UA 



E'a 



<Eb-E's) 



Hb 



(88) 
(89) 



The weak-interaction assumption is required for ((88)) to 
hold. Using ([55]) and and defining the heat and 
matter nonequilibrium constraints 



Ah = -Pa + Pb 
Am = Papa ~ PbPb 



(90) 



we find that (|87p can be expressed exclusively in terms 
of measured quantities Ea and ua 



R[E'A,n'A]EAjnA] ~ -Ah{E'A - Ea) - Amin-A - tia). 



Using jMl) and dSSl), we find 

S" = - > ' Pxi^rHxPa - pxT^rNxPo) 



E 

X=A,B 



(91) 



(92) 



S = - Pxi'TrHxm- Px'TrNxm)- (93) 



X=A,B 



From ([6T|) . the ensemble average of ([91]) is the time- 
integrated entropy production which has the familiar 
force- flux form of nonequilibrium thermodynamics (ssl . 

[TTolfi2i 



(94) 



(i?) « -A(TrH^p(t)-TriJ^p(0)) 
-^„(Tr7V^p(t)-TriV^/5(0)). 

The detailed FT follows from ^ and (l9T|l 

p{AEA,AnA) 



-{AhAEa+A^Aua) 



p{-AEA,-AnA) 



(95) 



Positive Ah [Am] means that Ta > Tb [PaPa > PbPb] 
so that the probability for an energy transfer AEa [of 
a particle transfer Aua] from A to is exponentially 
more likely than from B to A. 

Such a FT for heat has been derived in Ref. [9^. A 
simil ar FT for exchange of bosons has been derived in 
|104| |. This FT for pa rticles can also be derived from the 
OF of Ref. [Ulliel. Derivations of th is detai led FT for 
specific models are presented in section lIVB4l and lVB 21 



C. Steady-state fluctuation theorems 

We give simple qualitative and general arguments to 
show that the FT ([SS]) can be used to obtain a quantum 
steady-state FT for heat and matter exchange between 
two reservoirs through an embedded system. 
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We consider two reservoirs A and B (with Hamilto- 
nians Ha and Hb) each initially at equilibrium with 
its own temperature and chemical potential. A heat 
and matter exchange occurs between the two reservoirs 
through a weakly coupled embedded system (e.g. a 
molecule or a quantum dot). The total Hamiltonian is 
Htot = Ha+Hb+V, where V = Hs+Vas+Vbs contains 
the free Hamiltonian of the system Hs and the coupling 
term between each of the reservoirs and the system Vas 
and Vbs- The total Hilbert space is Ha x Hb x Hs- We 
use the index ix to distinguish between eigenstates of 
Hx with same energy Ex and number of particles nx, 
where X = A,B, S. We define the abbreviated notation 
i = {iA,iB,is) and a = {EA,nA, Eb^ub, Es,ns). The 
energy Ea and the number of particles is measured in 
reservoirs A at time zero and again at time t. We assume 

Po = Po Pa (Pa, Pa)Pb (Pb , P.b)Ps' {f3s , P's), (96) 

where p^' is the equilibrium system reduced density ma- 
trix. Since 

P[i' ,a';i,a] = \{i' ,a'\Ut\i,a)\'^{i,a\po\i,a) (97) 
P*^[^,a;^^a'] = ,a'\Ut\ha)\^i' ,a'\po\i' , a') (98) 



Eq. dnni) reads 



-(3A{iEA-E'A)-fiA{nA-n'A)) (99) 
-f3B{{EB - E'b) - MnB - n'g)) 
-Ps{iEs-E's)-^isins-n's)) . 



Since the system-reservoir couplings are weak, conser- 
vation laws of the total unperturbed system (Htot with 
Vas + Vbs = 0) implies that 

Eb-E'b ~ -{Ea-E'j,)~{Es~E's) (100) 
nB — n'g = n'^ — HA + n'g — ns - (101) 

This means that is equal to 

R[E'j^,n'^;EA,nA] w -AhiE'^ - Ea) - Amin'A ~ nA) 
+0{E's - Es) + 0{n's - ns)- (102) 

Since A and B are assumed macroscopic (i.e. reservoirs), 
the change in energy E'^ — Ea and matter ua — n'^ in 
reservoir A is not bounded. However, because system S is 
assumed small and finite, E'g — Es and n'g — ns are always 
bounded and finite. This means that in the long time 
limit, these contribution to R will become negligible in 
(UnSl). For long times, the FT ^SE^ with pI7^ becomes a 
universal (independent of system quantities) steady-state 
FT for the heat and matter currents 

1 p{AEa,ANa) a t ^ a t ^^r^',\ 

lim - In — — = Ahh + Arnirn , (103) 

t^oo t p[—AEa,—ANa) 

where = A£'^/t and /,„ = ANA/t are the heat and 
matter current between the system and the reservoir A. 
The r.h.s. of (|103p can thus be interpreted as an entropy 



production. A rigorous proof of (jlOSp has been recently 
given in Ref. [Sj]. In the long time limit, the steady-state 
FT (fTU5)) is similar to the detailed FT dM]). We note 
that the long time limit is related to the existence of a 
large deviation function (see appendix [C|) . We also note 
that when the system 5* is not finite, 0{E'g — £^5) and 
0{n'g — ns) terms in (|102p may not be negligible in the 
long time limit, as observed in Ref. JJW, 200]. Similar 
problems are expected if A and B are not "good" reser- 
voirs. A "good" reservoirs should allow the system to 
reach a steady-state. Since it is known that such reser- 
voirs cannot be properly described within the Hamilto- 
nian formalism, it should be no surprise that more sys- 
tematic derivations of quantum steady-state FT (|103p re- 
quire to use some effective (and irreversible) description 
of the embedded system dynamics. A common way to 
do this is the quantum master equation approach which 
consists in deriving an approximate equation of motion 
for the system reduced density matrix containing the ef- 
fects of reservoir through its correlation functions. As 
required for a "true" reservoir, the back-action of the 
system on the reservoir is neglected (Born approxima- 
tion). Such a derivation of the steady-state FTs will be 
presented in section ITVl [see (|14ip and (|149p ]. Another 
approach, is based on a system Greens functions descrip- 
tion. Here, the effect of the reservoirs appear through the 
self-energies. These derivations will be presented in sec- 
tion |VC] It has been recently suggested that finite ther- 
mostats (commonly used to model thermostatted classi- 
cal dynamics) could also be used to describe thermostat- 
ted quantum dynamics [Tlj . 



IV. HEAT AND MATTER TRANSFER 
STATISTICS IN WEAKLY-COUPLED OPEN 
SYSTEMS 

We now consider a small quantum system weakly inter- 
acting with a reservoir. Heat and matter exchanges are 
measured by a projective measurement in the reservoir. 
We will derive a generalized quantum master equation 
(GQME) for the GF associated to the system density 
matrix conditional to a given transfer with the reservoir. 
The statistics is therefore obtained from the solution of 
the GQME. When summing the GQME over all possible 
transfer processes, one recovers the standard quantum 
master equation (QME). 



A. Generalized quantum master equation 

We consider a single reservoir, but the extension to 
multiple reservoirs is straightforward. The total Hamil- 
tonian is the sum of the system S Hamiltonian, Hs, the 
reservoir R Hamiltonian, Hr, and the weak interaction 
between the two, V . 



H = Hq + V^Hs + Hr + V . 



(104) 
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We use the index s (r) to label the eigenstates of the 
Hamiltonian of system S (i?). The reservoir is initially 

assumed to be at equilibrium p^' — e^Pi^R^t^^R) j'E.ji. 

The measured observable is the energy Hn and num- 
ber of particle Nr in the reservoir. Since the measured 
observables commutes with the initial density matrix 
Po = ps{0)Pr, using we get 



G(A,i)=Trp(A,t), 



where A = {A,,, A^}, 



(105) 



(106) 



and 



= Ho + Vx. 



Obviously, p{t) = p{X = 0,t). 
We define the system GF 



ps(A,t) = TTRp{X,t) 



(108) 



which is an operator in the system space. Since psit) = 
PS (A = 0, t) is the reduced density matrix of the system, 
ps{X,t) is a reduced density matrix of the system condi- 
tional to a certain energy and matter transfer between S 
and B. Wc can now rewrite (|105p as 



G(A,i) = Tr5Ps(A,i) 



(109) 



We will derive a closed evolution equation for ps{X,t) 
by using projection operator technique and second order 
perturbation theory in V on p(A, t). By solving this equa- 
tion one can get G{X,t). Details are given in appendix 
Id] The final result reads 



/5s(A,t) 



{ 



-L[Hs,ps{X,t)] + ^Yl J (110) 

-Tr;^{VTT//(-T)p^V5(A,0} 
-Trfl{p^Vs(A,t)Vl^(-r)Vl«;} 
+Tr«{VTp^V5(A,t)T>;(-r)} 
+TTn{K{-r)p''^ps{X,t)V!^',}} , 



where 



e" 



(111) 



1. Generalized reservoir correlation functions 

We now consider an interaction of the form 

y^^S'^i?'", (112) 



where S"^ (-R") is a coupling operator of system S {B). 
It follows from (fT07)l that Va = where 

j^K ^ g^(Ahi?fi + A„7V/;)^Kg-^(AhiJ"ii + A„7VR) (113) 

For such interaction, (|110p becomes 

'ps{\t) = -l.[Hs,psi\t)] + ^Yl f (114) 

{-a«,.(r)^''^'''(-r)/55(A,t) 

-a,,,(-T)/55(A,i)5''(-T)5''' 
+a,.«(-A,-T)5''ps(A,t)^"'(-T) 
+a«.«(-A,T)^«(-T)p5(A,i)^«' } . 

Here we have defined the generalized reservoir correlation 
functions 



a««,(A,t) = TiRr'iR^^{t)R'^ 



(115) 



rr' ^^r' r 



where R'^{t) — e^^'^'^R'le"^^"* . The reservoir corre- 
lation functions are given by a^^'it) = ctnu'i^ = 0, i). 
For A = 0, (|114p therefore reduces to the non-Markovian 
Redfield QME of Ref. [U. 

The ordinary reservoir correlation functions satisf y the 
standard Kubo-Martin-Schwingcr (KMS) condition jll2| 



In the frequency domain 

dt 



(116) 



q;„k'(A,w) 



2tt 



e'-*a«„-(A,i) (117) 



the KMS relation reads 



(118) 



The generalized reservoir correlation functions satisfy the 
symmetry 

OiKK'H^h, Am}, t) = 

a^'^ii-Xh - i/3, -A„, + i/3/x}, -t) .(119) 

We note also that if i?" and S"^ are Hermitian, we further 
have 



a«K'(A,i) = q;*,„(-A, -i) 



(120) 



2. The Markovian and the rotating wave approximation 

Two approximations commonly used to simplify the 
QME may also be used on the GQAIE. The Markovian 
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approximation consist of setting the upper bound of the 
time integral in l|114p to infinity. The rotating wave ap- 
proximation fR WAI [2^ . [tHI (also known as se c ular ap- 
proximation [3^, Il56l | or Davis procedure [l03', 'l66l) is 
often used to impose a Lindblad form [23l Tl23. .166] to 
the Markovian QME generator in order to guaranty the 
complete positivity of the subsystem density matrix time 
evolution. Without RWA, the Markovian QME genera- 
tor can lead to a positivity breakdown for certain set of 
initial conditions due to small errors introduced on the 
initial short -time dyriamics by the Markovian approxima- 
tion [33,[M[iM[i0§,lii3,[i2l|- One has to note however 
that the use of the RWA is not always physi cally justi- 
fied and might miss important effects [30, Isil, 11091 Il71j . 
The RWA is equivalent to define a coarse-grained time 
derivative of the system density matrix on time s long 
compared to the free system evolution [3^ Il56l |. One 
easy way to perform the RWA consist in time averag- 
ing limy^oo 2T I-T '^^ generator of the QME in the 
interaction picture and in the system eigenbasis, using 



t:S{lu) ± iP— — lim 



1 



r)^0+ T] =P lUJ 



(121) 



Using these two approximation on the GQME (|114p . we 
find that coherences, Pss'{t) = {s\ps{t)\s') with s ^ s' , 
follow the dynamics 



where the relaxation rates are given by 



r , — 



- 27r(3KK'(0)S'^^S's 



(122) 



(123) 



and the modified system frequencies are 



OJ + UJs 



LU + LUss' 



.(124) 



The coherences evolve independently from the popula- 
tions [diagonal elements pss{i)\ and also independently 
from of each other. They simply undergo an exponen- 
tially damped oscillations which are independent of A. 
Populations, on the other hand, evolve according to the 
equation 



Pss{\t) 



(125) 



'^2T:a^^,{~ujss)S^,S^,Pss{\t) 
+27ra«,„(A, uj-ss)S';-^S';'^pss{\ t)} . 

The population dynamics depends on A. 



B. Applications to particle counting statistics 



We now calculate the particle statistics for different 
models and derive various steady-state FTs using the 
GQME. 



1. Fermion transport 

We consider a many electron quantum system attached 
to two metal leads which act as particle reservoirs. We 
shall denote the singe-particle eigenstates of the system 
and leads by indices s and i, respectively. The total 
Hamiltonian is H = Ha + Hb + Hs + V, where 



H 



X 



E 



(126) 



The coupling between the lead X — A,B and the system 

is Vx ^ Jx + Jx where Jx = Y.s,iex J^i^l^i and J^-f are 
the coupling elements between the system and the leads 
X. The total coupling is then 



V ^ Ja + Jb + Jl 



J 



B- 



(127) 



There is no direct coupling between the two leads, and an 
electron transfer is only possible by charging or discharg- 
ing the quantum system. The operators c{c^) represent 
the annihilation (creation) operators which satisfy the 
Fermi anticommutation relations 



clcl, +cl,cl 



0. 



(128) 



To connect with the notation of the Hamiltonian (|104p . 
we have Hji = Ha + Hb and V = Va + Vb- Apart 
from the difference in chemical potentials pA and pB with 
eV — PA ~ Pb, the two leads are assumed be identical. 

To count the change in the number of electrons in the 
lead A, the projection is done on A. Therefore (|llip for 
this model reads 



Vb 



e--^JA 



Vb 



(129) 



To get the second line, we used the relation Ja^a = 
{Na + 1)Ja- Substituting Eq. in (fTTU)) . the GQME 
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becomes 



i 1 /■* 

Ps{X,t)=--[Hs,psi\t)] + -^Y.J '^^ 

ss' ^ 

" {e^^at{-T)+al,{-r)}c,ps{XM{-r) 

+ {e^^ai.{T)+al.{T)}cA-r)ps{\t)c\ 
+ {e-'^/3,^,,(-r) +/3£,(-r)}ctps(A,t)c,,(-r) 
+ {e-'^(3t{T)+(}fArM{~r)ps{\,t)cs. 

- ass'{T)clcs'{~T)ps{X,t) 

- (3ss'{T)Cs'cl{-T)psiX,t) 

- ass'{~T)ps{X,t)cl{~T)Cs' 

- (3,s'{-T)ps{X,t)Cs'{-T)cl 



(130) 



where 



(131) 



are the equihbrium correlation functions for leads X and 
where a,,, (r) = a^, (r) + ag, (r) and (r) = (3f^, (r) + 

For A = 0, Eq. (fT30l) reduces to the QME derived in 
Ref. [88]. After applying the Markovian approximation 
described in section IIV A 21 (the upper limit of the time 
integral in Eq. (jl30p is extended to infinity), wc perform 
the RWA approximation which is equivalent to assume 
that the lead correlation functions are diagonal in s [ssl ] . 
Eq. (|130p then becomes 



Ps{X,t) = --[Hs,ps{X,t)] (132) 
+ E[ {e-''ai + ag}clPsiX,t)c^ 

S 

+ {e'^Pi+Pfjc.sPsi\t)cl 
- assCsclps{X,t) - PssclcsPsiX,t) 

The rates and are calculated by assuming a con- 
stant density of states a for the leads over the energy 
range around the Fermi level 



a 



X 



(3. 



X 



2tt 

2n 
If 



a\Jf\'il-fx{e.)) 
a\Jf\'fxies), 



(133) 



where /x(e) = [1 + e ^'^'^ ''^^] ^ is the Fermi function of 
lead X, and /3 = l/fc^T. These rates satisfy the relation 



^Ssr-'SS 



(134) 



The solution of (|132p allows to compute the time- 
dependent electron statistics between lead A and the 



system at any time. For A = 0, (|132p is the Lindblad 
QME derived in [88|] . Equation (|132p was first derived in 
Ref. [H^l by unraveling this QME. This means that the 
QME is interpreted as resulting from a co ntinuous posi- 
tive operator- valued measurement (23l . ll4l| on the system 
by the leads. This allows to construct probabilities for 
histories of electron transfers, and to use them to derive 
equations of motion for the GF associated with the prob- 
ability distribution of a net transfer of electrons during a 
given time interval, which are identical to (|132p . We thus 
find that the two-point projection method and the posi- 
tive operator- valued measurement lead to the same elec- 
tron statistics result in the weak coupling regime (with 
Markovian and RWA). A similar conclusion was reached 
in Rcfs. [46', 149]. 

In (|132p , the GF factorizes in terms of single orbital GF 
of the system, ps{X,t) — YlgLi Ps{X,t), where M is the 
total number of orbital and ps{X,t) is the single orbital 
GF, so that 



--e,[clcs,ps 



(135) 



{e-'^ai+afjclps{X,t)c, 
+ {e'V.1+/3f.R/5s(A,t)ct 

- UssCsclpsiX^t) - PssclcsPsiX,t) 



As discussed in Sec. (|IVA2p . the GQME p32)) . when 
expressed in the eigenbasis of the system describes an 
independent dynamics for coherences and populations. 
The coherences simply decay in time following damped 
oscillations while populations follow a classical rate equa- 
tion. If the eigenstates of each orbital are denoted by 
|ns), where = 0,1, the vector made of the pop- 
ulation of ps{X,t) in this basis denoted by ps{X,t) = 
{{0\ps{X, t)|0), (l|ps(A, t)\l)} evolves according to 



i^s{X,t)^T,{X)ps{X,t) 
where rs(A) is a 2 x 2 matrix 



r.(A)-( e-iV.+af. -Pss 



The eigenvalues of this matrix are given by 

ass + Ps 



ls±{X) 



(136) 



(137) 



(138) 



where 



/(A) = [e'^pf, + (3g){e-^^ai + afj + -(a.. - (Bssf. 

Since G(A,t) = llsGs{X,t), where G,(A,t) = 
{0\ps{X,t)\0) + {l\ps{X,t)\l), the long time limit of the 
cumulant GF is given by the dominant eigenvalue 

5(A) = hm ilnG(A,t) = V73+(A). (139) 

t— i-OO t ^ — ^ 
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Using p34l) and (fT38)) . we find that 7s±(A) 
7s-i-(— i/3ey — A), which imphcs that 



5(A) ^S{-i(3eV - X). 



(140) 



In appendix [Cl we show that this symmetry imphes the 
steady-state fluctuation-theorem 



Una = e''^^^ 

i^oo p[-~k, t) 



(141) 



where p(fc, t) is the probabihty of transferring a net num- 
ber k of electrons in time t from lead A to t he systern . 
Similar FTs have been derived in Refs. 0, [H, [isllTsot. 



2. Boson transport 

We consider a single oscillator mode at frequency eo//i 
Hs — ^oa\a(j coupled to two baths X = A,B at dif- 
ferent temperatures (3^^ and (3]^^ (fcs = 1) that consist 
in a collection of noninteracting bosons (e.g. phonons) 



Hr — Ha + Hb, where Hx 



The 



coupling is taken of the form V = Va + Vb, where 
Vx = J2iex Jio (^0 + a\){a\ + cLi). The subscript de- 
notes the system oscillator and i is for the i'th oscillator 
in the bath. J:^ is the coupling between the system and 
the i'th bath oscillator from X. All operators satisfy the 
boson commutation relations 



flso;, - a'g,as = dss', 
alal, - a'.al = agUs' 



0. 



(142) 



The system eigenstates have an energy iVgeo where Ns = 
1,2, - • •. We are interested in the statistics of the en- 
ergy transfers between the system and the A reservoir, 
so that the two energy measurements arc performed on 
system A. It can be shown that performing the RWA 
on the GQME is equivalent to assume from the begin- 
ning that the coupling term is of the simplified form 
= J2iex («1"o + a\ai). We thus have 



Vx = e^^^^ {Va + Vb 
= Ja(A) + j1(A) + 1/j 



e 2 



where 



Ja{\) 



(143) 



(144) 



ps{X,t) = --[Hs,ps{\t)] 



Substituting Eq. (fT43| in Eq. (fTTO)) , we get 

(145) 

-adaoalps{\,t) - aualdops{X,t) 

+ iate'^'° +a^)aopsi\t)^l 
+ {aie~'^'° + a^)alps{X,t)ao, 

where the rates au and ad correspond to the "up" and 
"down" jumps between the system states 

27rcr, 

2TTa 
- -^1-0 

1]~^ is the Bose distribution function 



|Jo^|2(l + nx(eo)) 



iJn^Pnx(eo). 



(146) 



"x(eo) 



and 



ad(u) 



''d(u) 



a^i N . The rates satisfy, 



d{u) 
B 



(147) 



For A = 0, (|145p is the Lindblad form QME derived in 
[8i,[i6l|. In the system eigenbasis {\Ns)}, Eq. (fT45l) de- 
scribes a populations dynamics which follows the equa- 
tion 



PNsi\t) 



- {ad{Ns - 



a^){Ns + l)pNs+i{\t) 
^l) + auNs}pNs{\t) 



y^) NsPNs-ii\t)XUS) 

where pxsi^t) = {Ns\psiX,t)\Ns) . Like Eq. ^ , 
(|148p may also be recast into a matrix form. However, 
unlike fermions, in this case since the matrix is infinite. 
p is an infinite dimensional vector and r(A) is a tridiag- 
onal infinite dimensional matrix. The determinant of a 
tridiagonal matrix can be expressed as a sum of terms 
where the nondiagonal terms always appear in pair with 
its symmetric nondiagonal term with respect to the di- 
agonal. With the help of Eq. p47p . this pair is sym- 
metric with respect to A — > — ieo(/3A — Pb) — A, so that 
dct{r(A)} = det{r(-ieo(/3A ~ Pb) - A)}. This implies 
that the eigenvalues have the same symmetry and there- 
fore that the following steady-state FT hold 

p{k,t) 



lim 

t^oo p[—k, t) 



eo{PA~PB)k 



(149) 



p(fc, t) is the probability that a net number of bosons are 
transferred from the reservoir A to the system i n a t ime t. 
Similar FTs have been derived in Refs. jSl,!!^!!!!. The 
transport statistics of bosons and fermions is different 
and was compared in Ref. f89'|. However, both satisfy 
the same type of FT [ (fTiTjl and (fT49ll ]. 



We have used ciiHA = (ci + HA)ai. Note that unlike 
fermions, Eq. (jl29p . in this case we have a factor in 
the exponential in the coupling, because we now mea- 
sure energy. However, in the present model the energy 
change is directly proportional to particle change, i.e. 
their statistics is the same. 



3. Modulated-tunnelmg 

In the above, fermion and bosons are transferred from 
one lead to another by charging or discharging an em- 
bedded system. We now consider electron tunneling be- 
tween two coupled leads, where the tunneling elements 
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are modulated by the state of an embedded system. Con- 
trary to the model of section IIVB 1[ the system never 
gets charged, however it affects the electron tunneling 
between the leads. This can happen for example if an 
impurity at the leads interface interacts with the spin 
of the tunneling electrons. The effect of this interaction 
is to modulate the tunneling elements between the two 
leads. This model of electron transfer was proposed in 
Ref. [148]. Here, we treat this model using the GQME 
approach. 

The Hamiltonian of the junction is of the form (|104p . 

where Hs is the system Hamiltonian and = Ha + Hb 

with Hx = J2iex ^iA^i = -S) ^'^^ the two leads 
Hamiltonian. The coupling between the two leads is of 

the form y = J + J'^, where J = X^ieA jes ■^ij^l'^j- The 
tunneling elements between the leads jjj = Jji are now 
operators in the system space. We measure the number 
of particles in the lead A. We then have 

V3, = e*^*^l/e-5^JV^ ^ ei^J + e-i^.P. (150) 

Substituting this in Eq. (IllOp . we obtain 

/}s{X,t) = -^[Hs,psi\t)]^ J2 (151) 

i€A,j£B 

' /^(e.)(l - /B(e,))4{4 W}ps(A,t) + h.c. 
+ fB{e,){l- fA{e.))Jl{Mt)}ps{Kt)+h.c. 

- /i3(e,)(l - /A(e.))e'' (4ps(A, t){4 (i)} + h.c) . 

- fA{e^){l - fB{e,))e-'^ ({^ W}ps(A, t)4 + /i.e.)" 
where 

{kit)} = ^J^ dre'^-^e-'^-Vye'^-^ (152) 

For A = 0, Eq. psip reduces to a Redfield equation for 
the reduced density-matrix of the system. A QME for 
the charge specific r educ ed density-matrix of the system 
was derived in Ref. [l48j . (|15ip is the evolution equation 
for the GF associated to it. 

When applying the Markovian approximation and the 
RWA to (|15ip in the system eigenbasis {|s)}, the popula- 
tions pss(A,t) = {s\ps{X,t)\s) evolve independently from 
the exponentially damped coherences according to 

Pss{X,t) (153) 

= {r,,,{x)p,,,,ix,t) - T,,,ix = o)pU\t)). 

s' 

The rates are given by 

r,,. (A) = e-'V,., + e'^e'^^^-'-^^^a,,- (154) 

where 

ass' (155) 

= ]^ E/^(^^)(l - fB{^,))MJ^^■'^')?5{^^, - Ess') ■ 



They satisfy the symmetry 

Ts's{-\ - i/3ey) = e'5^='^r,,,(A). (156) 

We define r(A) as the matrix generating the dynamics 
(|153p . Using Leibniz formula, the determinant reads 

N 

det{r(A)} = ^sgn(a)[]r..(.)(A), (157) 

a s 

where N is the order of matrix F and the sum is computed 
over all permutations a of the numbers {1, 2, TV}. 
sgn(cr) denotes the sign of the permutation, sgn(cr) — +1 
if a is an even permutation and sgn((T) = —1 if it is odd. 
Using Eq. (|156p , it can be shown that 

N 

dct{r(A)} = ^sgn((T)ne''''="<='r<,(,),(-A-i/3ey) 

a s—1 

N 

= ^sgn((T)nr.(.)s(-A-i/3eU) 

a s—1 

N 

= 5^sgn(a)nr..(s)(-A-i/3eT/) 

cr s—1 

= det{r(-A-i/3eU)}. (158) 

In going from first to second line, we used the fact that 
Yl^^^ ePE^.^s, ^ 1 due to Ef^, = 0. This prop- 

erty follows from the bijective nature of permutations 
which implies that for a given Esa{s) in the sum such 
that cr(s) — s', there will always be a in the sum 

that cancels the Eg'. Since the eigenvalues of F(A) sat- 
isfy the same symmetry property as the determinant, we 
get the same steady-state FT as (|14ip . where p{k,t) is 
the probability for a net number k of electron transfer 
from the lead A to the lead B. This shows that the FT 
(I14ip is not model-specific but rather a generic property 
of nonequilibrium distribution of electron transfers be- 
tween two leads. 



4-. Direct-tunneling limit 

When the system is decoupled from the junction, the 
tunneling elements between the two leads are given by 
Jij = Jij^- Using the Markov approximation, t — > cx) in 
Eq. (fT5^ . we get 

{4} = ^ ('^^(e. - e,) - iP^;^) (159) 

where is the principal part of x which we shall neglect. 
Under these approximations, it is possible to obtain the 
explicit form of the GF for the particle transfer statistics 
between the two leads. Substituting Eq. p59p in (|15ip 
and tracing over system degrees of freedom [Eq. (|109p ]. 
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we obtain 



(160) 



{fA{e^) + fsie,) - fA{e^)fB iej)}{cosX - 1) 

+ i{/B(eO-/A(ej)}sinA]G(A,t). 

The solution of this equation with the initial condition 
G(A,0) = 1 is 

G{X,t) = exp {tfiiie'^ - 1) + ti^i2{e-'^ - 1)} , (161) 

where 

2n 



fJ-i 



2n 



J2\J,^\^S{e.-e^)fB{e,){l-fA{e^)} 



We show in appendix [El that the probability distribu- 
tion associated to the GF (|16ip is a bidirectional Poisson 
process: the difference of two Poisson processes with mo- 
ments fii and fi2- Since the moments fii and fi2 satisfy 
/ii — e~^^^ fi2, the GF has the symmetry [see appendix 

El 



G{X,t) ^G{-X~iPeV,t). 
This immediately implies the FT 



Pi-k,t) 



^l3eVk 



(163) 



(164) 



which is satisfied at all times (transient FT) unlike (|141[) 
which only hold at long times (steady-state FT). The 
entire distribution p{k, t) is calculated in appendix [El 



V. MANY-BODY APPROACH TO PARTICLE 
COUNTING STATISTICS 



In previous sections, we formulated the counting statis- 
tics using a kinetic equation approach. This simple and 
intuitive approach makes some key assumptions. It as- 
sumes an initially factorized density matrix of the in- 
teracting systems so that initial Fock space coherences 
are ignored. Moreover, the approach is valid only in 
the weak coupling limit and it is not obvious how to in- 
clude many-body interactions such as electron-electron 
and electron-phonon. In this section we present a for- 
mulation of counting statistics based on superoperator 
non-equilibrium Green's functions (SNGF) [9l| which al- 
lows to relax these approximations. 



A. Liouville space formulation of particle counting 
statistics 

We consider particle transfer between two coupled sys- 
tems A and B described by the Hamiltonian 



H = Ha+Hb + V, 



(165) 



where the coupling reads V — J + . By choosing suit- 
able form for J, we can recover the different models stud- 
ied in section (jlVBp . For the present discussion, we do 
not need to specify the explicit form of J. 

The measurement of the net number of particles trans- 
ferred from j4 to -B is performed using a two-point mea- 
surement as described in Sec. [lT| Here the measured ob- 
servable is the number of particles in A. A measurement 
is done at time t = O.If right before this measurement 
the system is described by a density matrix |p(0):^, the 
measurement destroys all Fock space coherences and im- 
mediately after the measurement the density-matrix be- 
comes diagonal in the Fock basis. A second measurement 
is performed at time t. A difference of the two measure- 
ments gives the net number of particles transferred be- 
tween A and B. However if the particle transfer between 
A and B occurs though an embedded system, the two- 
point measurement of particle numbers in A measures 
the net particle transfer between A and the embedded 
system rather than between A and B. 

It will be convenient to work with superoperators in Li- 
ouville space [i,[6l[iS[9ll,[i33. These are defined in Ap- 
pendix[F| We shall denote Liouville space-superoperators 
by a breve and Hilbert space operators by a hat. Ha, Va 
and Hoa, where a — L,R, are the left and right super- 
operators corresponding to H, V and Hq — Ha + Hb- 
The probability of the net transfer of k electrons from A 
to B during the time interval t is [see Eq. (|G6p ] 

p{k,t) = «^|4-fet^(i,0)P„|p(0)», (166) 



where U{t,0) — e is the time evolution operator 

in Liouville space and _P„ is the projection operator asso- 
ciated with the measurement of n electrons in A. |p(0)3> 
is the interacting density matrix when the counting starts 
and contains coherences in the number operator basis. It 
is constructed by switching on the interaction V from the 
infinite past, where the density-matrix |p(— oo)^ is given 
by a direct product of the density-matrices of systems A 
and B, to t — 0. 



|p(0)>= C//(0,-cx))|p(-oo)>, 



(167) 



where 



Ui{0, — oo) = exp^ |--y drV2F_(T)| (168) 

with V2V1(t) = Vi(T) - Vr{t) [see Eq. ([F8])] and 

VaiT) = ul,{T,0)VMT,0), a = L,R (169) 
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where 

;7o(t,0) = 9{T)e-r^^"o~^. (170) 
The GF associated to p{k, t) is defined by 

G{X,t)=J2e'^''pik,t). (171) 



We define 



Substituting Eq. (fTTTj) in Eq. p66)) . we get (see Ap- 
pendix [G| 

G{\t)= -^G(A,A,i) (172) 



(2^)^ 



with 



G(A,A,i) 



|p(-oo)>, 



(173) 



where V^h{t)) = ^^(tlW) - Vfl(7flW) with 7L(t) = 
e{t){K + A/2) and 7K(t) = 0{t){K - A/2). The GF (fT7511 
includes the initial < = correlations between systems 
A and B in the density matrix. These correlations are 
built through the switching of the coupling V from t — 
— oo and i = 0. In the absence of such correlations, the 
initial density matrix is diagonal in the number basis and 
G(A,i) = G(A,A = 0,i) [i.e. p(0) commutes with Na]- 
Below, we show how G(A, A, t) can be computed. 



B. Electron counting statistics for direct-tunneling 
between two systems 



We next apply Eq. p73p to calculate the electron cur- 
rent statistics for the direct tunneling model of section 
pVB4p . The Hamiltonian is given by (|165p . where 



J — Jijclcj, 



(174) 



where J*j = Jji. Hamiltonian Ha aud Hb are general 
and can include many-body interactions. The exact form 
for Ha and Hg is not necessary in the present discus- 
sion. A noninteracting electron model, as studied in Sec. 
pVB4p . will be considered in the next subsection. 

We now define the superoperators J, and TV corre- 
sponding to the operators J, and the number operator 
Na for the system A. These satisfy commutation rela- 
tions 



[Jl,Nl] = -Jl , [JINl] = jI 



(175) 



Using these commutation relations in Eq. (|G14p , we can 
write 

V^a(7a(0) = exp{-i7„(t)}J„ -|-exp{i7„(i)}4 . (176) 



Z(A,A,t) = lnG(A,A,i). 



(177) 



Expanding the time-ordered exponential in (jl73p wc can 
compute the GF and the cumulant GF perturbatively in 
the couphng Jab- Since <^I\V-\p{—oo):s>= 0, to second 
order we obtain 



Z(A,A,t)==-^ / / UT2 



dTi / dr9 

— oo J —oo 



(178) 



«/|ry_(7(Ti),ri)Vl(7(T2),T2)|p(-(»)» . 
Substituting Eq. P75)) in pTg]) we get 

Z(A,A,<) =Z(o)(A,t)+z(i)(A,A,t) (179) 

where 

(A, A, t) = 2(e'^/2 - e-'^/2)Re{e'^M^(i)} (180) 

(181) 

and 

Z^"\\,t) = (c-'^ - l)W^^lit) + (c'^ - l)T4^i°i(0 (182) 

are the contributions coming from time evolution from 
t = — oo to i — and from t = to time i, respectively, 



and W{t) EE W^\{t) - W\'^{t) with 
^Asit) = ^J^ dhj^ dt2«/|4(ii)Jl(t2)|p(-oo)» 

W^BAit) = '^^l / dt2 ^I\JL{tl)jUt2M-^> 

^isW = J^J «/|Jfl(tl)jl(t2)|p(-00)» 



.(1), 



^BAit) ^ '^^l y ^^^2 «/| Jl(<i)4(^2)|p(-«^)» . 

(183) 



1 



From P77|) and pT^l) we get 

G(A,A,t)=e^*'"(^^*)e2*"(^'^^*). (184) 
Substituting this in (|172p . the GF is obtained as 

G(A,t) = G(")(A,t)G(i)(A,t) (185) 

where 

G("Ha,0 = exp{z(°)(A,t)} (186) 

G^'Hx,t) ^ £''^exp{Z«(A,A,t)} (187) 

The cumulant GF is finally obtained as 

Z(A,t)=Z(")(A,i)-)-ln / — exp{Z«(A,A,t)}. 

(188) 
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The second term on the rhs of Eq. (|188p is the contribu- 
tion due to the initial correlations that exist between sys- 
tems A and B right before the first measurement. When 
these initial correlations are ignored, i.e. initial density 
matrix is a direct product of the density matrix of A and 
B (or equivalently [iV^, /5(0)]=0), Z^^) = 0. 



1. Effects of initial correlations 

Here we discuss the corrections to the electron statis- 
tics due to correlations between A and B in the initial 
density matrix. We show that these contributions do 
not affect the first moment (the current) but only higher 
moments. 

Using (jlSOp and expanding in A, we find that 



exp{Z«(A,A,i)} 



Tl = 



J2 ^ sin" ( 9 ) [e-'^'Wit) + e^^W*(t)]" 



E 



n(>fe),fc=0 



(2»)" . „ fx 



xW'''''{t)W''*{t)e-'^^"~^''\ (189) 
Integrating over A, (|187p becomes 



GW(A,i) = 1 + E ^imt)P"sin2" 

n=l ^ 



(190) 



By differentiating (|185p with respect to A, we can fac- 
torize the moments in two parts, (A:")o which does not 
depend on the initial correlations and A*^") which does: 



where 



(fc"(i)) = (fc"(i))o + A(")(i), 



Qn 



(191) 



dX 

n OQ 21 



A=0 

2l-m+k( 1 ^fc 2/ 



k=l 1=1 m=0 



' ( -i\k ZLf-i n/-i 

x(/fc"-'=(t))o|T^(t)P'. (192) 



"Cfc = n\/{k.\{n — k)\) are the binomial coefficients. 

We find that H^S^^t) — 0, i.e. initial correlations do 
not contribute to first moment, which is the net current 
from A B. However, they do contribute to higher 
moments. The correction to the second moment is 



A(2)(i) = -32|T4^(t)p 



(193) 



We see that initial correlations always tend to decrease 
the second moment. 



2. The thermodynamic limit 

We consider now the limit where A and B can be as- 
sumed to have continuous spectra. We treat them as 
non-interacting electron leads and show that initial cor- 
relations do not contribute to the long time statistics. 

This corresponds to the model discussed in section 
IIV B 41 In this limit, the rates Wab and Wba given in 
Eq. p83p can be calculated explicitly. The Hamiltonian 
for two systems {X — A, B) is 



(194) 



Using the fact that the density-matrix at < = — cxd is a 
direct product |p(— oo)^= ®\p^b^^ 



</|JK(ri)j|(r2)|p(-oo)> 



(195) 



,-(-ri-T2) 



«/|JL(ri)4(r2)|p(-c^)» (196) 
- E \J^J?fB{^j){^ - /A(Q))e"^-(^^-^^)} , 

where ujij = e — Cj and fx{() = (exp {/3(e - ^x)} + 1)""^ 
is the Fermi function for the system A{B) with /i^ and 
/is denoting the chemical potential of systems A and B. 
Remembering that 

Jo \ ^ij/^ 



and that 



dri I aT2 
la 



(77+ T iwjj)^ 



using (|182p we find that 

Z(°) (A, t) = (e-'^ - l)fi2t + (c-^ - l)Aiit , (197) 

where fii and fi2 are given by ()162p . G'"'(A) is therefore 
identical to the GF for a bidirectional Poisson process 
obtained in (fTeTIl within the GQME. 

The rate W{t) which appears in the expression for Z*^^^ 
in Eq. (fT^ reads 



W{t) = ^|J,, 



(/a(wj) - fsiuj,))- 



-i{uJij-iri+)t 



(198) 



Taking the continuous limit of the leads' density of states, 
we fi nd th at for long times W{t) becomes time indepen- 
dent [Tel]. Therefore 

lim -Z{X,t) = lim -Z^"\X,t) , (199) 

t — ^00 t i — ^00 t 

which shows that the long time statistics is not affected 
by the initial correlations between A and B. 



5(A) 
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Electron counting statistics for transport 
through a quantum junction 



We next apply Eq. (I172|) to calculate the current 
statistics in the transport model of section H V B 1 1 where 
a quantum system (e.g. a molecule, chain of atoms or 
quantum dots) is embedded between two much larger 
systems A and B. Notice that here the two-point mea- 
surement of the particle number in A does not measure 
the net particle transfer between A and B as stated in 
section IV Al but rather the net particle transfer between 
A and the embedded system. The particle transfer statis- 
tics for this model was studied in section lTVB ll using the 
GQME approach. Here, we express the transfer statistics 
in terms of the SNGF ^90,|9l| of the quantum system. By 
connecting this powerful many-body formalism with the 
two-point measurement, we can study more complicated 
models. The effect of eigenbasis coherences in the quan- 
tum system (which requires to go beyond the RWA in the 
GQME approach) and the effect of many-body interac- 
tions in the quantum system can be easily incorporated 
into the SNGF approach via the self-energy matrix. In 
presence of many-body interactions, the SNGF theory in- 
volves a self-consistent calculation for the Green's func- 
tions together with their self-energies. This goes beyond 
the weak coupling limit of the GQME. The simple form 
for the lead-system interactions (|127p allows us to obtain 
analytical results for the corresponding self-energy and 
hence the GF. Electron-electron interactions will provide 
an extra (additive) self-energy matrix computed in Ref. 

The Hamiltonian of the model is given by (|126p and 
P27p . The superoperators Hoa and Va corresponding to 
Hq = Ha + Hb + Hs and Va(B) can be obtained by using 
Eqs. (|FIT|) in Eqs. and p?7)) . We get 



(^xRC'^R 



(200) 



xeA,B,S 



and 



where 



Vo, = JA,a + JAa + Jb,c. + JI^ (201) 



Jx,L = ^ Jsicli^CsL 

s.iex 

Jx,R = J^^^sRc\ji. (202) 

s,i£X 

The superoperators Jx,l and jj^ ^ satisfy the commuta- 
tion relations [9l| 

[Ja,l,Nl\ = -Jal; [Ja,r,Nb] = Jam. (203) 



[I 



A,L'- 



= J 



A,L' 



[J 



A,R 



(204) 



and [Na, Jb,. 
we obtain 



0. Using these in (|G14p . 

(205) 



+ exp{i7„(t)}J^,„(i)-H Jb 



J 



't 



where Jx.a = Jx,a{la — 0). Note that in (|205p . expo- 
nential factors are associated only with superoperators of 
the lead A. This is because the measurement (projection) 
is done only on A. 

We can now use (I205P in p73p to compute the GF. 
|p(— oo)3> in (|173p is given by the direct product of equi- 
librium density-matrices of the system and the leads, 



|p(-oo)> = \ps))®\pa))®\pb)) 



\Px)) = 



» 



(206) 
(207) 



where and S^; are respectively the chemical potential 
and the partition function for system x. 

Using Grassmann variables and a path-integral formu- 
lation, the GF (|172p can be expressed in terms of the 
Green's functions of the quantum system. In Appendix iHl 
we present a derivation in terms of Liouville space super- 
operators. For a Hilbert space derivation see Ref. [l02| |. 
Some useful properties of Grassmann variables used in 
the derivation are summarized in Appendix H] The final 
result for the GF, Eq. pTg)) . is 



G(A,A,i) 



(208) 



with 



Z(A,A,i) = 



/ rfrlnDct [g-\T = 0) - S(r,r,7(r))] , 

J — OO 

(209) 



where g{t — t') and E(i,f') are Green's function and 
self-energy (due to system-lead interaction) matrices in 
v,u' = — representation. The Green's function matrix 
satisfies 



gZ {t-t')^d{t-t')6,,,5ss' (210) 



and the self-energy matrix is 

srr;(t,t',7(i),7(i')) = 



E 

X ii'ex 



where 



4+ (7) = J-s- 

for i ^ A and 

^++(7) 



(7) = JU^- 
(7) = Me'-'- 







-e'''«)/2 
e'''")/2 (212) 



(213) 
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for i E B. One important point to note is that while 
(t, t') (zero-order system Green's function without in- 
teractions with leads) is causal and g ^{t, t') = o(90l.[9ll|. 
this is no longer the case for and E ^ which depend 
on 7. This is due to the fact that when 7^ ^ 7/j, the ket 
and the bra evolve with a different Hamiltonian. The 
cumulant GF is then given by 



Z{\t) 



In / _G(A,A,i). 







27r 



(214) 



Equation ([2T4|l with ((208ll and (|209|) give the statistics for 
the net particle transfer between lead A and the quantum 
system embedded between A and B. 



1. Long-time statistics 

At steady-state all the two-time functions, such as 
g{t, t') and S(i, i'), depend only on the difference of their 
time arguments. We factorize time integration in Eq. 
(j209p in two regions, one from —00 to and other from 
to t. Since 7(t) = for negative times, Eq. (jGlSp . we 
obtain 



Z(A,A,t) = Go 



(215) 



[ drhiDet [.g~i(r = 0) - E(r, r, 7(t))] 



The term Gq , which is independent on time and 7 comes 
from integration < = — ootoi = and contains all initial 
correlations between system and the leads. Substituting 
for the self-energy (|21ip . we notice that since the matrix 
elements and J^^ appear at the same time, the A 
dependence drops out. We can recast (|209p for long times 
as 

^ InDet [g-^{Lo) - Y.{uj, A)] + Go.(216) 

At long times the first term in (|216p dominates, and the 
current GF is given solely by the first term in (|216p . 

5(A) = hm -Z{\t) 

t^oo t 

= \nDct[g-\oj)-j:{uj,\)] 

^ 1 g lnDet[x-i(u;)]. (217) 

Thus, as in Sec. IVB 2| we can conclude that contribu- 
tions coming from the initial correlations between the 
system and the leads do not effect the long-time statis- 
tics. 

We shall compute the self energy in frequency domain. 
Since the leads are made of non- interacting electrons, 
their zeroth order Green's functions in frequency domain 
are 

9^ri.) ^ .i^(.)H.-]lH 

gT-+{uj) = -2T:iSu'{2Muj)~l)S{hw-e,). (218) 



Substituting this in (|21ip . the self-energy matrix in the 
wide-band approximation is obtained as 



S++(a,,A) = 



2 



(219) 



E-7(c.,A) = --rf,, 



(220) 



Si7(o;,A) 



-^Ti,[{e''-l)il-U{u;) (221) 



^ -iTfA2fBiu^)-i) + ^rt 

X [(e'Vl)(l-/A(c.)) 

-(e-'V 1)/A(a;)] , 



(222) 



where F^, = '2TrJ2iex JisJs'iS{LO - e,). 

Note that when A = 0, ~ as it should 

be (causality) [iOl , and E , E++ and E ^ reduce to 
usual retarded, advanced and correlation (Keldysh) self- 
energies, respectively. 



S-7(a;) 



--F , 



E++(c.) = -F, 



E-+(^) - iF,,, - 2{ri,fA{co) + Ff,,/s(^)) (223) 

where F = F^ -HF^. 

The retarded Green's functions for the molecule is then 
given by 



1 -1 



(huj — e)i — i- 



(224) 



where 1 is the identity matrix and = [R ]^ is the 

advanced Green's function. 

Finally, we transform the self-energy matrix from the 
— (E) to the L, R (E) representation. This can be 

achieved by the matrix transformation, E — Q^^Y,Q [90l |. 

where 



^/2 I 1 -1 



This gives the matrix EQ,^(a;, A) with elements 



(225) 



i{Ti,Uto)+rgjB{co)) (226) 



^Tl{^) = -^r^s' + i{ri'fA{cj) + rgjBicj)) (227) 

Ef;,(c., A) = iFf,,/B(c.) + iri,fAiu^)e-'^ (228) 
E||;(c., A) = -iFf,,(/B(c.) - 1) - iTfAfAico) - l)e^\ 

(229) 

where the A dependence occurs only in Ei/j and E^^. 
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Equation ((2T7|) together with ((2T8)) and ((2T9l) - (|222)) 
gives the long-time current statistics within the two-point 
measurement approach. It contains the full information 
about the coherences in the system eigenbasis through 
the self-energy matrix E and can therefore be used to 
study effects of coherences on the current statistics. 



2. Recovering the generalized quantum master equation 



The GF (piT)) is different from the GF obtained using 
the GQME approach (|139[) . We are now going to show 
in what limit ((2T7l) reduces to (fT39l) . 

Assuming that the S matrix is diagonal, Y,ss' = 
Sss'^-ss, the determinant = \g~^ - S] in Eq. (PTT]) 

factorizes into a product of determinants correspond- 
ing to each orbital s, |x~^l = lis \x7s\^ and Z{X,t) — 
J2s ^s{X,t) becomes the sum of GF for individual or- 
bitals. We note that the assumption of a diagonal E 
matrix amounts to ignoring the coherences in the quan- 
tum system eigenbasis and is therefore the analog of the 
RWA in the GQME approach. In the following we com- 
pute Zg. For clarity, we omit the orbital index s in the 
self-energies. Since from ([ ^ - ((^ ^ll = [^rr]*, we 
can write 



\x7l\ = {u;-esr-\f:LL\' 



LR^RL- 



(230) 



Substituting this in (|217p . we get for the long time 
cumulant GF 



^^(A)^ / ^In iu;-esr + \^LL\''-^LR.^RL 
Zn 



(231) 



In order to compute the frequency integral we first obtain 
the A-dependent current by taking the derivative with 
respect to A 



duj 



dxi^LR^^RL) 



27r (w-e,)2 + |SLi|2-I]iflI]«i 
Using (1226]) - ([229]) ■ we get 

/,(A) = ir^r^ 

dLO }b{u;)[1 - /A(w)]e'^ - ./a(c^)[1 - IbH] 



(232) 



(233) 



27r 



M{\lj) 



where 
M{\,uj) 



r^ + r-^r^ [/b(c.)[i - /A(t^)](e'^ - 1) 



+ /A(c.)[l-/s(c.)](e-'^-l)] 



(234) 



Assuming that the couplings with the leads are weak 
ksT » so that resulting broadening is small com- 
pared to Eg, the contribution to the integral comes mainly 
from the center of the Lorentzian. This allows us to 
replace w = es in the Fermi functions inside the inte- 
grand. We therefore need to consider the poles lu = 



Ss ± iy/ M(A, e^). Computing the residues at the poles, 
we get 



^.(A) = 



-iF^r^ 



2v/M(A,e,) 



fB{es)[l - fAicsW^ (235) 
-/A(e.)[l-/B(e.)]e-'^ 



Since Zs{X) = -^(A'), we finally get that 



^.(A) = - 2 + 



(236) 



which coincides with the GF obtained from the GQME, 
(fT39l) . 



3. The Levitov-Lesovik formula 



Equation ([^17)) with (PT5| - ((^^ is the most general 
formula for the transport statistics at long times for a 
system of non-interacting electrons. It includes the ef- 
fects of coherences between the various tunneling chan- 
nels (system orbitals) available to an electron tunneling 
between the two lead. This is due to the non-diagonal 
structure of the self-energy in the Hilbert space of the sys- 
tem, Eqs . (I219 [)-(| 222D . Here, we recover Levitov-Lesovik 
formula [l2Ct Il2ll | for the counting statistics. For that 
we again assume diagonal self-energies. As discussed in 
the previous subsection (|V C 2[) . the cumulant GF in this 
case is simply the product of the GFs for each orbital. 
Thus all orbitals contribute independently to the electron 
transport. 

Using self-energy expressions (|226p - (|229l) . the GF 
(|23ip can be expressed as 



5. (A) 



dio ^ 
2^ 



-In 



hr^r^[/B(c.)(/A(^)- 

h/A(c^)(/B(^)-l)(l 



(237) 



Using (|224p . we can write for orbital s 

Substituting (l^gS)) in (^57)) . we obtain 

5s(A) = -2 J ^\n\R--iu;)\ 
duj 



(238) 



+ J ^l^{l + ^H[/BM(/A(^)-l)(l-e'^) 
+ /a(^)(/b(c.)-1)(1-c-^)]} (239) 

where T{lj) — r^T-^\R'^g^ {uj)\'^ is the transmission co- 
efficient for the tunneling region. The first term on the 
r.h.s. of (|239p can be ignored since it does not contribute 
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to the average current or its fluctuations (independent on 
A). Therefore 

Ss{x) = y'^hi{i + rH[/BH(/^H-i)(i-c'^) 

+ /^H(/BH-l)(l-e-'^)]} (240) 

which is the Levitov-Lesovik formula [l20l Il2ll . Il22| . It 

has been recently generalized to a multi-te rmin al model 
for a non- interacting tight-binding model |l58l |. Equa- 
tion (|240p is valid to all orders of the coupling. The only 
approximation required to obtain the Levitov-Lesovik ex- 
pression p40p is to ignore the coherence effects between 
different orbitals in the tunneling junction. Notice that if 
T(a>) is small, we can expand the logarithm in Eq. (|240p . 
This is equivalent to making a perturbation in the cou- 
pling V. The leading order in the expansion gives (|197p 
with 11621) . 

Since fA{uj)[l - /b(^)] = e'5^^/s(c^)[l - /a(w)], it is 
straightforward to see that the GF (|240p satisfy, 5(A) = 
S{~\ - iPeV), and the FT ([141]) follows. 

Taking the derivative with respect to A of the GF (j240|) 
at A = 0, the average current is 



that a net amount of energy or matter k crossed the sys- 
tem during a given time. The cumulant GF defined as 



then possesses the symmetry 

Z{X,A) ^ Z{iA- X,A) 



(242) 



(243) 



Taking the derivative with respect to A of both sides and 
using (|242p , we find that in the ^ ^ limit 

— [Z(A, 0) - Z(-A, 0)] = -i^2(A, 0). (244) 
The cumulant GF is expressed in terms of cumulants as 

(iA)^ 



Z(A,^)= ^i^if„(^). 



(245) 



m— 1 



Using (|245|) in (I244p . we find at each order in A, that 

[1 - (-1)"] ^^™(0) = Km+m. (246) 



duj 
2^ 



(241) 



which is the Landauer-Buttiker expression for the average 
current through a tunneling junction with transmission 
coefficient T{uj) (l6j . 



VI. NONLINEAR COEFFICIENTS 

As we have seen, the FT implies a specific symme- 
try of the GF which depends on the nonequilibrium con- 
straints imposed on the system. For weak constraints, 
i.e. close to equilibrium, this symmetry can be used to 
derive fluctuation-dissipation re latio n as well as Onsager 
symmetry relations [§, [sl, [H, Ill7l |. A systematic ex- 
pansion of the GF in the nonequilibrium constrains al- 
lows to derive similar fundamental relations further away 
from equilibrium. This has been done for stoc h astic 
systems for counting statistics [H, [HI [HO, [Hi 
and for the work FT [6]. FTs therefore provide a sys- 
tematic approach for studying generalized fluctuation- 
dissipatio n re l ation s such as previously considered in 
Refs. m Eli 1113 • 



A. Single nonequilibrium constraint 



We assume that a FT of the form p{k,A) — 
e^'^p{—k^A) holds in a system maintained in a nonequi- 
librium steady-state by a single nonequilibrium con- 
straint A, where p{k,A) is the probability distribution 



Equation (|246p implies that at equilibrium, odd cumu- 
lants are zero and event cumulant are related to the 
derivative with respect to the nonequilibrium constraints 
of the nonequilibrium odd cumulants when approaching 
equilibrium 

^2m-l(0) = (247) 
if2™(0) = 2^/^2^-1(0). (248) 

Below we show that this leads to the well known fluctu- 
ation dissipation relations. 

We next consider the second derivative with respect to 
A of both sides of Using and and after 

some algebra, we find in the ^ ^ limit that 



dA^ 



[Z(A,0)-Z(-A,0)] 

6)2 



(249) 



Using (|245p . we find at each order in A that 



[i-(-ir 



dA^ 



i^m(O) 



(250) 



= [i + (-i)™+i] ^/W(o). 



This relation is only useful for odd m and implies 

-^K2^-m=QjK2.n{Q). (251) 

This procedure can be continued for higher derivative of 
Z{iA — X,A) with respect to A. 
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We can always expand the average process in term of 
the nonequilibrium constrain as 

Ki{A) = Ki{0) + l'^^'^A + L^^^A'^ + 0{A^). (252) 

i^^) is the Onsager coefficient. Using (|247l) . ([248]) and 
(pgTj) for m = 1, wc find that Ki{0) = and that 



^(2) 



_d_ 

OA 
1 



2dA^ 



ifi(O) 



K2{0) 
2 

- lA 
~ 2dA 



K2i0). 



(253) 
(254) 



([253]) is a fluctuation-dissipation relation. As an illus- 
tration, we consider a biased quantum junction such 
as in section IIVB II k represents the number of elec- 
tron crossing the junction and the nonequilibrium con- 
straint is given by ^ = PeV, where V is the potential 
bias across the junction. In this case, close to equi- 
librium, (/) = /Se^yL^^^ is the average electrical cur- 
rent through the junction and 6^/^2(0) is the Fourier 
transform of the equilibrium current correlation func- 
tions at zero frequency. (|253p indicates that the resis- 
tance of the junction, which characterize the dissipation, 
is related to the current fluctuation at equilibrium by 
R = dv{I) = /3e2i(i) = /3e2i^2(0)/2. 



B. Multiple nonequilibrium constraints 

When multiple nonequilibrium constrains are applied 
to the system, the FT can be used to find important 
symmetries of the response coefficients 0, Q • In case of 
N nonequilibrium constraints, the cumulant GF reads 



Z({A^},M^}) =ln J2 ^''M{k,},{A,}) I ,(255) 
\{k-,} 

where X-k = k^^-y- We assume that it satisfies the 

FT symmetry 

Z({A^}, {Ay}) = Z{{iAy - A-,}, {A^}) . (256) 



Proceeding as in section IVI Al we find that (|244p gener- 
alizes to 



d 



dAfi 



[Z({A,},{0})-Z({-A,},{0})] 



(257) 



d_ 

'dXp' 



-i— Z({A^},{0}) 



The cumulant GF can be expressed as 
Z{{Xy},{Ay}) 



(258) 



{m^} = l ^7=1 ^' 



where the cumulants read 



zm,{A^}) 



(259) 



The generalisation of ((2461) is found using ([258| in ([2571) . 
so that at a given order in the A's 



N 



1 - ]^(-l)"'^ 

7=1 



^^{n.-,}({0}) 

dAp 



(260) 



If we choose {rriy} = {S^a}, we get that 



dAp 



K{5,^+5,,}m) ■ (261) 



Close to equilibrium, the average processes can be ex- 
panded in term of the nonequilibrium constraints as 



({-^7}) — ^ ^ La^A^ ~\~ ^ ^ L(x'y~f' A^A^ 

7 7>7' 



(262) 



Since the (Onsager) linear response coefficients are given 
by 



i, 



Q/3 



dAp 



(263) 



using (|261[) . we find the Onsager reciprocity relation 

L^p = L/3a . (264) 



The generalisation of (|250p to multiple nonequilibrium 
constraints reads 



52 



dAadAi3 
92 



[Z({A,},{0})-Z({-A^},{0})] 



(265) 



-Z({A^},{0}) 



92 



_dKdAp~'' " 
This implies that 



dXpdAa 



(266) 



7 

dK 



({0}) 



dAc. 

N 



1 _Y[ (^\>n-,+S,,\^^{]Jh±b^}!^^ 

^ ) dAp 



7=1 



For {m^} — {5~f0}, we get 
5'^{5..}({0}) 



Leap ~ 



dAadAfi 



(267) 



,(dK{s,,+s,,}m) , aif{,^,+5,„}({0}) 



dAa dAfj 
which implies the expected symmetry Lgap = ie/3Q 
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VII. CONCLUSIONS AND PERSPECTIVES 

The approach to quantum statistics adopted in this 
review is based on a two-point projective measurement. 
This, together with considerations about the symmetry 
between the forward and the time-reversed quantum 
dynamics, allow to recover from a simple and unified 
perspective all previously derived fluctuation theorems 
(FTs) for quantum systems (transient as well as steady- 
state FTs). This was the object of section HIl and Hill 

A generalized quantum master equation (GQME) is 
presented in section IIVI for a quantum system weakly 
coupled to reservoirs. It describes the evolution of the 
generating function (GF) associated with the system 
density matrix conditional to the outcome from a two- 
point measurement (of energy or number of particles) 
on the reservoir. When summed over all the possible 
outcomes, the quantum master equation (QME) for 
the system reduced density matrix is recovered. This 
formalism has been applied to various model systems and 
used to directly demonstrate the validity of steady-state 
FTs. 

The GQME formalism circumvents the unraveling of 
the QME, used to calculate the quantum statistics of 
particles or ei iergy, and origiiially d eveloped in quantum 
optics [H m, Halli, [m \lM EMI ■ Since the unraveling 
of a QME is not unique, a continuous time measurement 
on the reservoir is assumed in order to connect the 
resulting quantum trajectories to measurable quantities. 
This procedure is only possible for Markovian QME 
which preserve complete positivity [in the rotating wave 
approximation (RWA)]. In this regime, the GQME 
formalism predicts the same statistics as the unraveling 
formalism. This equivalence between the two types 
of measurements in th e weak coupling limit was first 
found in Ref. [i^, Il49j | . This results from the fact that 
the reservoirs are assumed to always remain described 
by the same canonical or grand canonical equilibrium 
density matrix (sil . [s^ and are therefore not affected by 
the measurement. The net number of particles or the 
net amount of energy transferred during a given time 
interval is then the same if the reservoir is continuously 
monitored or only measured twice at the beginning and 
at the end. The unraveling of non-Markovian QME has 
been an active fie ld of research during this last decade 
psi . [sol . Il69l . ll7Cll |. but the connection between the re- 
sulting quantum trajectories and measu rable quan tities 
is not straightforward [H, 113, [zl, M, [US 02^- In 
the GQME formalism, the connection to measurable 
quantities in the non-Markovian regime is unambiguous. 
Exploring non-Markovian effects on the particle or en- 
ergy statistics could be an important future application. 

In order to go beyond the approximations used in 
the GQME formalism (i.e. initially factorized density 
matrix, weak coupling), we presented an alternative 



approach based on superoperator non-equilibrium 
Green's functions (SNGF) in section |V| This Liouville 
space formalism provides a powerful tool for calculating 
the particle statistics in many body quantum systems. 
Using this formalism, we showed that initial coherences 
in the basis of the measured observable do not affect the 
steady-state counting statistics and the FT. This is to 
be expected since at steady-state, the long time limit 
destroys the information about the initial condition. We 
showed it using a non-interacting electron model for 
both direct and indirect (transport) tunneling between 
two reservoirs. However, for transient FTs such as the 
Crooks relation, the assumption that the system density 
matrix is initially diagonal in the basis of the measured 
observable seems unavoidable for the FT to be satis- 
fied. We applied the SNGF formalism to compute the 
counting statistics in some simple models and discussed 
the limit in which the statistics predicted by the QME 
is recovered. The Levitov-Lesovik formula for electron 
tunneling between two reservoirs, which goes beyond 
the weak coupling limit of the QME, was also recovered. 
We discussed the approximations required to recover 
the Levitov-Lesovik expression from a more general 
result expressed in terms of the SNGF for the tunneling 
region. In particular, we showed that when several 
energy channels are available to tunneling electrons, the 
Levitov-Lesovik approach does not capture the quantum 
coherence between different channels. This amounts to 
ignoring the off-diagonal elements of the self-energy in 
the eigenbasis of the system. 

Transient FTs (valid for arbitrary time) have been 
presented in IIIIBI The work FT derived for isolated 
driven system in section IIII B II is always valid since, 
besides an initial canonical density matrix, no assump- 
tions have been made. The work FT for open driven 
system derived in section IIII B 21 assumes an initially 
factorized canonical density matrix between the system 
and the reservoir and a definition of work which is only 
consistent for weak system-reservoir interaction. The 
transient FT for direct heat and matter transfer between 
two finite systems and derived in section fill B 31 assumes 
that the systems are each initially at equilibrium and 
weakly interacting. The steady-state FTs (only valid for 
long time) presented in section IIII CI and derived more 
systematically in section IIVI assumes a weak system- 
reservoir coupling and the RWA. However, the FT has 
been recently shown (numerically) to hold for QME 
without RWA [ml and the Levitov-Lesovik formula 
presented in section [V C 31 is obtained nonperturbatively 
and satisfies the FT. FTs seem therefore to characterize 
universal feature of nonequilibrium fluctuations in 
quantum as well as in classical systems. 

We now discuss some future perspectives. 

We mentioned in the introduction and in section 
im that an alternative approach to counting statistics. 
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where the GF used is an influence functional following 
from a path integral description of the system-detector 
interaction, has been developed during the last decade. 
It is only in a semi-classical limit that the two-point 
measurement approach predicts the same statistics as 
this approach. Determining the region of applicability 
of both prescriptions is an open problem that could lead 
to a better understanding of quantum measurements. 

Various numerical methods have been developed for 
using the Jarzynski relation to efficiently calculate equi - 
librium free energies of classical systems |ll8l IllOl Il85l | . 
Extending these methods to quantum systems will be of 
interest. 

Finally, we note that in this review we have focused on 
systems maintained in a steady-state distribution by a 
single non-equilibrium constraint. Investigating systems 
subjected to multiple nonequilibrium constraints could 
reveal interesting features. 
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APPENDIX A: TIME-REVERSED EVOLUTION 

We explain why ([M)) corresponds to the time-reversed 
expression of the two-point probability ([¥]) and discuss 
how to physically implement a time-reversed evolution. 
The effect of a static magnetic field is also discussed. 

In order to implement the time-reversal operation in 
quantum mechanics, it is necessary to introduce the an- 
tilinear operator Q {Qi = —iQ) which satisfies 8^ = 1 
(i.e. = 6) [HO, [Hll- An arbitrary observable A 

can be even or odd with respect to the time-reversal op- 
eration, i.e 

eie = caA , (Al) 

where — ±1. For example, the position operator R 
is even {en — 1) while the momentum P or angular mo- 
mentum L are odd (ep.L = ^1)- It can be verified that 
the Heisenberg commutation relations are preserved un- 
der the time-reversal operation. When acting on a time 
dependent Hamiltonian H{t; B) that depends on a static 
magnetic field B, wc get 

eH{t;B)e ^ H(t;-B) . (A2) 



If a forward evolution operator [as in ([2]) but with a static 
magnetic field] evolves according to 

jU{t,Q-B)^-'-H[t-B)U(t,Q-B) , (A3) 

with the initial condition i/(0,0;i?) — 1, than the time- 
reversed evolution operator is defined by 

J7t,.(t,0;-B) = eU{T -t,0;B)U\T,0;B)e 

= eU{T ~t,T;B)e , (A4) 

and its evolution is given by 

^^Ut,{t, 0; -B) = -^H{T - t- -B)Utr{t, 0; -B) , (A5) 

with the initial condition C/tr(0, 0;i3) = 1. This can be 
verified using the change of variable t ^ T — t u\ (jA3[) . 
then multiplying the resulting equation by Q from the 
left and by [/^(T, 0;i?)O from the right and then using 
(IMl) and ((Xl . 

From now on we choose t = T (the time at which the 
time reversal operation is performed is t), and we define 

p{t) = U{t,0;B)poUHt,0;B) (A6) 
8/5*^(06 = UtAt,0; - B)e p'J eulit, 0; -B) .{A7) 

We note that by multiplying (|A7|1 by 8 from the left and 
from the right, we get 

p''^{t)^U\t,0;B)p'JU{t,0;B) . (A8) 

We verify that if p^' = p{t), then p*''(i) = po- This 
means that, as for classical systems, if a system initially 
described by po evolves according to the forward evolu- 
tion between and then the time-reversal operation is 
applied and the resulting density matrix is evolved ac- 
cording to the backward evolution during a time t and 
finally the time-reversal operation is again applied, the 
resulting density matrix is the initial condition /5o. It fol- 
lows from this discussion that if the two-point probability 
[with a static magnetic field B\ is defined as 

P[auao\= (A9) 
Tr [Pa, U{t, 0; B)Pa„poPa„ {t, 0; B)Pa, } , 

the time-reversed expression of this two-point probability 
has to be defined as 

P''[ao,at]= (AlO) 
Tr {PaMt, 0; -B)Pa,epllQPa,uUt, 0; -B)P,„} . 

We note that we could have included the final time- 
reversal operation in the definition, but it has no effect 
anyway due to the trace invariance. By inserting 8^ in 
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between all the operators in (|A10|) . and using (|A4|) with 
T = t, we find that 

P''[ao,at]= (All) 



which is identical to the definition used in (p4)) . It is con- 
venient to use (jAlip as a starting point because it allows 
to avoid mentioning the presence of a static magnetic 
Field. However, it is important to keep in mind that the 
physical evolution corresponding to the time-reversed dy- 
namics associated to a forward dynamics with an Hamil- 
tonian H(t; B) is an evolution with an Hamiltonian where 
the driving protocol is time-reversed, where the sign of 
the static magnetic field is changed H{T — t; —B) and 
where the initial condition is Qp^qQ. 



APPENDIX B: FLUCTUATION THEOREM FOR 
COARSE-GRAINED DYNAMICS 

Here, we show that using a coarse-graining of the ini- 
tial density matrices, R defined in section HH AI becomes 
a measurable quantity and (R) a difference of Gibbs - 
von Neumann entropy. We follows closely Refs. [28l . ll95l |. 

We define 

R[at,ao] = In = -R''[ao, at] (Bl) 



and 



p{R) = J2 P[at,ao]5iR~ R[at,ao]) 

at, do 

p"^{R) = P''[ao,at]S{R~R''[ao,at]) . (B2) 



Note that (|B1|) . in contrast to ([50)) . is expressed exclu- 
sively in terms of measurable quantities (eigenvalues of 
A{t)). An integral FT follows 

(e-^) = J2 «o]e-«['^-'^«l = 1 , (B3) 



which implies (i?) > 0, as well as a detailed FT 



(B4) 



The coarse-graining of a density matrix p within its non- 
measured part reads 



EPa A 



(B5) 



where Pa = TrpPa is the probability to measures a, and 
da is the number of states with the value a. When, as in 



15l|, such a procedure is applied to p^ and po, {R) can 
be related to an entropy change. In this case 

P[at,ao] = TT{U\t,0)Pa,Uit,0)Pa„}^ 
P^'iao.at] = Ti{U^{t,0)Pa,U{t,0)Pa„}^ . (B6) 



(B7) 



Therefore, using (jB6p in (|Bip . we get 

R[at,ao\ = s|j' - Sao , 

where 

tr 1 Pat 1 

The average of R now reads 

(i?) = J2 ao]P[at,ao] ^ S'' - S , (B9) 



Pao 



(B8) 



where 



(BIO) 



are the Gibbs-von Neumann entropies associated to the 
coarse-grained density matrix p^ and po. Indeed, if the 
coarse-grained density matrix p is used in the expression 
for the Gibbs-von Neumann entropy S = Trplnp, we get 

S = Y^a^aPa- 



APPENDIX C: LARGE DEVIATION AND 
FLUCTUATION THEOREM 

Below, we briefly describe large deviation theory and 
show that a symmetry of the long time limit of the 
cumulant GF such as (|140p or (|158p translates into a 
steady-state FT for the probabilities. 

We consider a probability distribution p{t,k), where 
/c is a counting variable associated to a continuous time 
random walk (we assume that the waiting time distribu- 
tions have a finite first and second moment). For fixed 
time, the central limit theorem is only valid up to a given 
accuracy in a central region of the probability distribu- 
tion hows width does not converge uniformly with time. 
Large deviation goes beyond the central limit theorem 
and allows t o describ e the behaviour of the tail of the 
distribution [l65l Il8l| . It relies on the assumption that 
the probability p{t,£,) that ^ — k/t takes a value in the 
interval + d£] behaves as 



p(t,0 = C(f,i)e^(«^' 



(CI) 



where the large deviation function (LDF) is defined by 



i?(a^ lim llnp{t,0 



(C2) 
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and where 



lim -\nC{£_,t) = . 



(C3) 



We will show that the LDF is determined by the long 
time limit of the cumulant GF given by 



5(A) = lim -lnG{t,X) , 

t^ca t 

where the moment GF is defined as 



G(t,A)=^p(i,fc)e 



(C4) 



(C5) 



Note that for convenience, we have absorbed a factor — i 
in the definition of A compared to the standard definition 
of the moment GF used in the main text. The GF can 
be rewritten in terms of p{t, ^) as 



We can then rewrite (jC6P as 



(C6) 



(C7) 



At long times, the main contribution to this integral 
comes from the value of ^, 5*, that maximizes the ar- 
gument of the exponential. £,* is therefore the value of 
^ such that A = ^|^^^.. At long times, using steepest 

descent integration, (|C7p becomes 
G{t,X) 



(C8) 

1 1 d^mo I ^ I (f f.2f 



eWr)-Ar)tc(r,i) 



de 



t 

2^ 



We assumed concave to have a maximum. Substi- 
tuting (1C8| in gives 



where 



5(A) = i?(C) - AC , 



A = 



(C9) 



(CIO) 



This shows that 5(A) is the inverse Legendre transform 
of the LDF. By taking the derivative of (|C9p with respect 
to A, we get 



dS{X) _ dR{O d^ d^ 
dX d£_ dX dX ^ 



which using (|C10[) leads to 



dSjX) 
dX 



(Gil) 



(G12) 



This shows that the LDF is given by the Legendre trans- 
form of 5(A) 



R{0 - 5(A) + Ae 



(G13) 



By taking the derivative of (IG12P with respect to A and 
using the derivative of (|C10|) with respect to ^, we can 
confirm that i?(^) is concave because 5(A) is convex. 

We now assume that the cumulant GF satisfies the 
symmetry 



5(A) =5(A- A) 



(G14) 



We note that the symmetry (|G14|) with the standard def- 
inition of the moment GF would read 5(A) = S{iA — A)]. 
Using the symmetry (IC14I) . Eq. (|C13|) implies that 
R{~0 = S{A - A) - (A - A)$, so that 



Using Eq. (|Gip . we get 



lnm = Aet + ln.^(^'^) 



Using (jCSp . this gives the steady-state FT 

hmiln^i^^AC, 
t-+oo t pit, —k) 

which is often written as 

p{t,k) t_»oo Ak 

— e 



p{t,-k) 

Eqs. (|14ip and (I149p are of this form. 



(G15) 



(G16) 



(G17) 



(G18) 



APPENDIX D: DERIVATION OF THE 
GENERALIZED QUANTUM MASTER 
EQUATION 

Eq. (|106p satisfies the equation of motion 

'p[\t) = L^p{X,t) = -'-{H^p[X,t) - p{X,t)H^>.) (Dl) 
= {U + vL\)p{X,t) 

= -'j[Ho,K\i)]'^\{Vxp{X,t)-p{X,t)V->) , 

where we multiplied U by a scalar v to keep track of the 
order in the perturbation expansion below. Superoper- 
ators are denoted by a breve [see appendix |F] . In the 
interaction representation where 



pi{X,t) 



-Lot-/ 



p{X,t) = e*"°'p(A,t)e-^^«* , (D2) 



L\{t) = e-^°*i',e^°* , 
Dip takes the simple form 

'pi{X,t)^vL'^{t)pi{X,t) 



(D3) 
(D4) 
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By integrating Eq. (jD4p and truncating it to order , 
we get the perturbative expansion 



pi{\,t) 



W{\t)p{Qi) = e-'^«*e^^*p(0 
Wo{\,t) + vWi{X,t) 

WW2{\t) + 0{v'^) p(0) 



(D5) 



where 



Jo Jo 
The inverse of W{t) reads 

W-\X,t) = Wo{X,t)-vWiiX,t) 



(D6) 



*dT r dTL'^{T)L\{T-T) 



(D7) 



Indeed, one can check that 1^(A, i)T^~i (A, t) = i+e'(i;-'^). 
For later use, we also notice that 



W{X,t)AW^^{X,t) = vWi{X,t)A 



(D8) 



W2{X,t)A~Wi{X,t)AWi{X,t) +0{v^). 



We define the projection superoperator (acting in reser- 
voir space) 



(D9) 



where is the equilibrium density matrix of the reser- 
voir. We used the Liouville space notation [see appendix 
|F]. P satisfies the usual properties of projection superop- 
erators P + Q ^ 1, P'^ = P,Q'^ ^ Q a.nd PQ = QP = 0. 
When acting on the density matrix p{t), the projection 
operator gives 

P|p(A,t)» = |p5(A,t)»®|p"«'» . (DIO) 

We now let P and Q act on the density matrix of the 
total system in the interaction picture (|D5P and find 

P|p/(A,t)» = PT^(t)(P + Q)|p/(0)» (Dll) 

Q|p/(A,t)» = QWit){P + Q)\pi{0):$> . (D12) 

Hereafter, we consider initial conditions such that 
Q\p{0)^= 0. This means that the reservoir part of the 
initial condition is diagonal in the reservoir eigenbasis 
and is thus invariant under the evolution when v = 0. 
Taking the time derivative of Eq. and Eq. (|D12p 

and using |p/(0)>= W^"^(A, t)|p/(A, i)>, we get 

P|p/(A,t)> = PW{X,t)PW-\X,t)P\pi{X,t)X^l3) 
+PW{X, t)PW-^X, t)Q\pj{X, i)> 

Q|p/(A,t)> = QW{X,t)PW-\X,t)P\pi{X,t)M^14) 
+QW{X,t)PW-'\X,t)Q\pi{X,t):$> . 



So far these equations are exact. If we restrict ourselves 
to second-order perturbation theory in v, we can obtain 
the important result that the P projected density matrix 
evolution is decoupled from the Q projected part. Indeed, 
with the help of Eq. (jD8|) , we have 



PW(X,t)PW-'^{X,t)Q = vPWi{X,t)PQ (D15) 
WPW2{X, t)PQ - v'^PWi{X, t)PWi (A, t)Q + 0{v^) . 

The first two terms of the right-hand side are zero be- 
cause PQ = and the third one also because 

P^i{X,t)P = J2 |PK'»«rr|L^(A,i)|p^'»«rV|(D16) 

vanishes since commutes with Hji. 

Having shown that the relevant projected density ma- 
trix evolves in an autonomous way, we will now evaluate 
the generator of its evolution using second-order pertur- 
bation theory. Again using Eq. (jD8p . we find that 



PW{X,t)PW-H\t)P = vPWiiX,t)P (D17) 
WPW2iX, t)P - v'^PWiiX, t)PWi{X, t)P + 0{v^) . 

The only term of right-hand side which is not zero is the 
second one [see Eq. (|D16p ] whereupon we get 

P|p/(A,i)» = v'-P I dTL'^{t)L'^{t ^ T)P\pi{X,t)-> 
Jo 

+0{v^) . (D18) 

Now leaving the interaction representation and using the 
fact that Pe-^°* = e~^^*P, we obtain 



P|p(A,t)» = L5P|p(A,t)» 



(D19) 



+w^e^^'P / drL^(<)L';^(t-T)e-^^'P|p(A,i)> 
By taking the trace of (|D19p we get 



PS 



(A,t) = Lsps{X,t)+v^y] f dT 

^ Jo 



(D20) 



X e^-* «rr|L',(t)iA(t-T)|p^^'»e-^-*p5(A,i) . 
Explicit evaluation leads to Eq. (jllOp . 



APPENDIX E: BIDIRECTIONAL POISSON 
STATISTICS 

The GF of section llV B 41 corresponds to a bidirectional 
Poisson process. We give here some basic properties of 
this process. 
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The GF of the probabiHty distribution p{k) can be 
expanded in terms of moments (fc") as 



(El) 



The Poisson distribution and its GF are given by 



p{k) 



fi e 



fc! 



G(A) =exp{^(e'-'-l)} . (E2) 



Note that /i ~ (fc). If fc = fci — fc2 where p(fci,fc2) = 
Pi(fci)pi(fc2) and pi(fc) and P2{k) are Poissonian, we get 
that 

G(A) = Gi(Ai=A)G2(A2 = -A) 

= cxp{Aii(c'^-l) + M2(c-'^-l)}. (E3) 

If the average of the positive process is related to the 
average of the negative one by = ^2 exp(— A), we find 
that the GF displays the FT symmetry G(A) = G{A-\). 
By inverting Eq. (jE3|l . we get 



p{k) 



-(Ml+M2)gAfc/2j 



(Mi - M2) 
sinh(yl/2) 



(E4) 



where Ik is the modified Bessel function of order fc. 

APPENDIX F: LIOUVILLE SPACE AND 
SUPEROPERATOR ALGEBRA 

In LiouviUe space, a N x N Hilbert space operators p 
is mapped into a vector |p» and a superoperator A 
(linear map) acting on an operator p becomes a N'^ x N'^ 
matrix acting on the vector \p:^: Ap ^ [ol. Isil [90l. 

[oil . l20l| . We recall some basic definitions 

scalar product : <A|B>= Tri^B , (Fl) 
identity : 1 = ^ \nn':$>^nn'\ , (F2) 

n.n' 

|rm'>^ \n){n'\ , <Cmi'| ^ \n'){n\ . (F3) 

Useful consequences of these definitions are 

<nn'|nn'> = SnnSn'h' (F4) 
<nn'|A> = {n\A\n) (F5) 
<1|A> = Tri (F6) 

We define left and right LiouviUe space operators as 

Al\x:^^Ax, Ar\x:^^ xA. (F7) 

We also define 

A+ = ^iAL + Au), A^ = ^{AL-An). (F8) 

This linear transformation is symmetric. The inverse 
transformation can be obtained by simply interchang- 
ing + and — with L and R, respectively. Thus most of 



the expressions in the following are symmetric and the 
indices used to represent superoperators can take both 
+, — and L, R values without any other change. The 
advantage of the +, — representation is that a single op- 
eration A_ in LiouviUe space represents the commuta- 
tion with A in Hilbert space. Thus all the intertwined 
commutations, that appear in perturbation expansions 
in Hilbert space transform to a compact notation that is 
more easy to int erpr et in terms of the double sided Fyn- 
mann diagrams |133l | . Similarly a single operation of A+ 
in LiouviUe space corresponds to an anticommutator in 
Hilbert space. 



yl_|X> ^ 



A+\X > 



1 

1 



{AX - XA) 
{AX + XA) 



(F9) 
(FIO) 



For any product of operators in Hilbert space, we can 
define corresponding superoperators in LiouviUe space 
using the following identities. 



{A^Aj---Ak)L = A^LAjL---AkL 

{A,Aj---Ak)R = AkR---AjRA,R. 

Applying this immediately gives, 
1 



(Fll) 



iA.A,)- 



{A^A,)+ = 



[■A^+,AJ+ 



2V2 
1 



] + [A-,A,^] 



(F12) 



2^/2 



{A+,Aj+} + {A,^,A,^} 



[A+,Aj_ 



] + [A^,A^+] 



(F13) 



Equations (|F11[) - (|F13P are useful for recasting functions 
of Hilbert space operators, such as Hamiltonian, in terms 
of the superoperators in LiouviUe space. 

Another useful c^antity in LiouviUe space is the time 
ordering operator T; when acting on a product of super- 
operators (each at different times) , it rearranges them in 
increasing order of time from right to left. 

fA,^{t)Aj^{t') = Aj^{t')A,a,{t), t<t' 

= i„(t)ijV3(i'), t' <t.{FU) 

where a, f3 — L,R,+,—. Note that, unlike the Hilbert 
space where we have two time ordering operators describ- 
ing the evolution in opposite (forward and backward) di- 
rections, a LiouviUe space operator T always acts to its 
right and therefore all processes are given in terms of for- 
ward times alone. This makes it easier to give physical 
interpretation to various algebraic expressions commonly 
obtained in perturbation expansions which can be con- 
verted readily in terms of different LiouviUe space dia- 
grams. 

We finally note that using (jF6P and (jFTJ we get for 
a = L, R that 



{A)=Tt{Ap} 



(F15) 
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and using jFTl, and (|F15|) . we get 

<C/|i_|p>=0 , < /|i+|p>= V2(i) . (F16) 



APPENDIX G: PROBABILITY DISTRIBUTION 
FOR ELECTRON TRANSFERS 



In the model considered in Sec. IV Al we consider 
electron transfer between system A and B. We mea- 
sure the number of electron in system A at time and 
time t. The number operator for system A is defined 
as iV J2i<£A^i^i' where c^(c) are creation (annihila- 
tion) operators. Only the coupling V can induce electron 
transfer: [Ha + Hb,N]=0 

The total density matrix follows a unitary dynamics in 
Liouville space 

|p(t)»= ;7(i,o)|p(o)»= c/L(t,o);7^(t,o)|p(o)» ,(Gi) 

where 

U{t,0) = exp{-iV2i7_t} (G2) 

with is the superoperator corresponding to the total 
Hamiltonian, \/2H^ = Hl — Hfj and 

Ua{t,0) = cxp{-iV2iJ„0, a = L,R. (G3) 

By measuring the number of electrons in A, when the sys- 
tem right before the measurement is described by |p(0)3>, 
we get the outcome n with a probabihty <C/|-Pn|p(0)^ 
and the density matrix of the system after the measure- 
ment becomes P„|p(0)3>, where the projection operator 
in Liouville space is defined as 

(G4) 



P. 



SK{n- NL)SK{n~ Nr) 
2^ dXdX' 



g-iA(n--Afi)g-iA'(n-ArR) 



5k is the Kronecker delta and Na are the left and right 
superoperators corresponding to the number operator in 
A. We have P„P„' = Sxin — n')P„ and 



exp{iAA^Q}P„ — exp{iAn}P„ . 



(G5) 



The net number of electrons k transferred between A 
and B during time i is a fluctuating quantity. The prob- 
ability for measuring k electrons during this time interval 
is given by 



p{k,t) = «^l^n-fcC/(i,0)P„|p(0)» 



(G6) 



Substituting (|GTj) and (|G4l) in (|G6l) and using ((G5|) 
with the fact that left and right superoperators commute, 
we get 



</|e'^i*^C/L(i,0)e-'^i*^ 

xe'^=*«[/jj(t,0)e-'^=^*«|p(0)> 



(G7) 



Making the change of variables, Ai = —A — A/2 and 
A2 = A — A/2, we get 



f rl\ 

where the GF reads 



G(A,i) 



dA 
2^ 



G(A,A,t) 



(G8) 



(G9) 



and 



G(A,A,t) = «/|e-'(^+^^/2)Jv^^^(^^o)gi(A+Ax/2)JV^ 

^gi(A-Ax/2)JV«jj^(^^0)e-'(A-^^/2)*«|p(0)» . (GIG) 

Equation (|G9p is identical to the trace of p{X, t) defined 
in ([25]). 

The density matrix right before the first measurement 
{t = can be constructed by switching the interaction V 
adiabatically from the remote past, t —00. This gives 

G(A,A,i) =«/|C/o(t,0)C//(7(i),t,-oo)|p(-c»)»(,Gll) 

where 

C//(7(t),t,-(X3) = exp+j-iy" dr\/2y_(7(r),r) 

Uo{t,0) = 9{t)exp{-iV2{Ho_t)}. (G12) 
We define 

V2V-{j{t),t) = Vi(7L(T),T) - Vfl,(7fl(r),r) (G13) 
where 

VLilL{T),T) = e-'-^^M^^ (Vl(t)) e'^^(")*^ 
Vfl(7fl(T),T) - e'^«(^'*« (vH(r)) e-'^«(-)*« (G14) 

with Va = Ja + J a ^ud 



7L(i) = ^(i)(A + A/2) 
jR{t) = e{t)(A-X/2). 



(G15) 



The time dependence of operators in (|G12[) is in the in- 
teraction picture with respect to Hq. 



Equation ((GTT|) is the GF used in Eq. P75)) . 



(G16) 



APPENDIX H: PATH-INTEGRAL EVALUATION 
OF THE GENERATING FUNCTION FOR 
FERMION TRANSPORT 

The fermion coherent states \ip) are defined through 
the eigenvalue equation for the Fermi destruction opera- 
tors, Cx\ipi) = ^xiW and (V'j|4 = (M^xi, where V and 



31 



ijj are independent Grassmann variables (see Appendix 
IJ) which satisfy antic ommutation relations similar to the 
Fermi operators [l40| |. 

It is convenient to introduce coherent states in Li- 
ouville space corresponding to the superoperator Cxa, 
X — a,b, s, as 



CxL\i')) = V'a^LlV')) 



(HI) 



The state j-!/")) can be expressed in terms of the vacuum 
state 



IV')) 



3SI^(-'/'xtC;^-l/JxRCa;Ji) I 



and 



IpE^CV'xiCxi-l/JxRclH) 



(H2) 



(H3) 



Note that is not the hermitian conjugate of cn [9l| . 
Grassmann variables ipa and anticommute between 
themselves and with the creation and annihilation oper- 
ators. Note that, unlike usual fermion case, we now have 
four generators for the Grassmann algebra, two corre- 
sponding to each index a. Using (jHip - (jH3p . it can be 
shown that 



((^'I^CflJV')) = 

((^'kWlV')) = 

((v^'icicLiv)) = 

{{^l>'\cnch4>)) = 



V'flV'L (('/'' I V")) 



(H4) 
(H5) 
(H6) 
(H7) 



These matrix elements will be useful in the path-integral 
formulation below. The scaler product of two coherent 
states is 



Grassmann variables satisfy the closure relation. 



1= / i?(v;^)c-2:.v;..„^...|^^))((^^l 



(H8) 



(H9) 



where V^ip-tp) = Il,,x.aidtpixa)(d'ipixa). 

We next switch to -I-, — notation[90j, represented by 
the index v — +,—. Superoperators Hq- and V_(7) are 



V2Ho- - 5]e.(4^ 



V2V^m,t) = ^^4^,<(7(t))cLWc..'((PlO) 

where J^^, (7)(= J^Jx) is 2 x 2 matrices for v,v' ~ 
with elements 

^j.t(7) = J-,7(7) = J..'(e^"^+e^^«)/2 
•/.V(7) = J-t(7) = J..'(e^^^-e*^-)/2 (Hll) 



while for x,x' = 6, s, J^x'^l) = ^xx'^l) = ^xx' and 
•/.V(7) = J.".t(7)=0. 

We will encounter the matrix element of an exponential 
operator of the type 



((V^'le 



m = ((^'le^-^^-lV-)) 

= ((^'|e2:x^-(^xL2xi,-2x«2L)|^)) 



where in going from second to the third line we used 
(jH6p and (|H7p . We can now make the linear transfor- 
mation from L/R variables to the +/— variables. In 
Hilb ert space this corresponds to the Keldysh rotation 
[l02j | . Using this transformation we can write above ma- 
trix element as 

((Vy|e^*°-|V)) = e^-'^'<"'''^-'^-'((i/;'|i/;))(H12) 

where = — tx and e^^ — e^^ = 0. This matrix 
element (jII12l) can also be obtained directly by formally 
defining the Grassmann variables corresponding to +, — 
operators, Cxi, and c|,^, by tpxv and 'ipxv, respectively, and 
using (jH10|) . We shall use the +/— formulation in the 
rest of the section. The advantage of using this notation 
is that we directly work with the retarded and advanced 
functions which are naturally linked to the observables 
(when A = 0). 

We can express the trace in Eq. (jGlip in terms of the 
coherent states basis, 

G{X,A,t) = 

(H13) 



We next divide the time from to t in Eq. (|II13|) into 
N equal segments of length 5t and introduce the closure 
relation (jH9p after each time interval. We then get. 



G(A,A,t) = 

v{H^){(^o\pi~^))){{i^i\Uih{t),t,-^)\^Po)) 

N 



(H14) 



Here the index i on ipi carries time index so that ipi+i is 
at St time ahead of ipi . 

JJi can be formally evaluated by dividing the time in- 
terval from the initial time — tg (at the end wc can put 
to — * oo ) to i in N' number of equal time steps. We then 
get 



((^i|C//(7(t),t,MlV'o)) 



(H15) 



i|e 



l^o)). 



Here > is small enough so that only the linear order 
term contributes. The exponential can then be factorized 
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into products of exponentials. By inserting the identity 
between exponentials, we obtain (repeated indices are 
summed over), 

((^,+i|e-'^^-(^^(*-)'*')*"*'|V4) (H16) 

The second matrix element of the evolution operator 
Uo{t) between two coherent states is 



'**'((V',+i|V'»)lH17) 



Using Eqs. pn6|) and pTTTl) in Eq. pTTS)) . we obtain 
for the GF, 

M=N+N' 

G{\K,t)^ v[i^mMp{-^))) n 



exp <^ iipix^ I t^ipixu I m 



SU 



exp 



(H18) 



Here -/(U) = for i < M - N. 

Setting M ^ oo, to ^ oo and SU 0, we get 



G(A,A,t) = / I?[V;V']e*'^('^''^) (H19) 



where in the continuous time notation Vlipip] = 
YIt action S{'ip,ip) is defined as 

- ^^ArW^i;: {l{T))i>x'u'{T)) . (H20) 



g^'^i is a 2 X 2 matrix corresponding io v,v' ^ which 
satisfies 



«^ - ) 9xx' (t, t') - Sit ' t')Sx,x'S,,' . (H21) 



Using the integral identity for independent Grassmann 
variables ?7, 77, S and n 

J 2?(ry,r;)e-''"^'^"^e''''''+*'''^' = Dei[A]e'''['^Ii'''^(H22) 

we can trace out the leads' degrees of freedom to obtain, 

Gi\,A,t) ^ J V[i;ilj]e'^'-'f"f'^ (H23) 

with 

si4>^) = J dTdT'^,ur)[g::r\r,T') 

-Sr;;(r,r',7)l ^.v'(r'). (H24) 



The self energy S(7) is 

Ks'(t,t\j)= (H25) 
E E (7W)5mi,i')^;^/(7(i')), 

xx'^A.B 

where repeated arguments are summed over and gxx' 
are the Greens functions for the non-interacting leads. 
The counting parameter appears in the self-energy only 
through couphng terms V^^ (7). Finally, using Eq. 
(|H22p we can perform the Gaussian integral in Eq. (|H23p 
to obtain, 



G(A,A,t) = exp[2:(A,A,t)] 



(H26) 



where 



Z(A,A,t)= / dTlnDet[5-i(r = 0)-S(r,r,7(r))]. 

J —00 

(H27) 

Here gir^r') and S(r, r',7(r)) are matrices in -|-,— su- 
peroperator indices and defined in the system space. This 
result for the GF was used in (|209p . 

APPENDIX I: GRASSMANN ALGEBRA 

Here we briefly review come properties of the Grass- 
mann algebra used in Appendix |h1 Fermion coherent 
state s 1 77) are defined in terms of the vacuum state |0) 

EM. 



{V\ = {0\e^'^ = (0| + (Oh*c 



(11) 
(12) 



where 77 and rj* are two independent complex numbers. 
Here we consider a single degree of freedom. This can 
be generalized easily for several degrees of freedom for 
which, |?7) = e^'^''''^\0). 

Since coherent states are the eigenstates of the annihi- 
lation operator, c\rj) — ri\rj), from Eq. (|IT|) . we have 



(13) 



which is a consequence of — (c*)^ — 0. Also since 
c, anticommute, it can be shown from the eigenvalue 
equations 



7777* + 77*77 = 0. 



(14) 



The independent variables 77 and 77* which satisfy Eqs. 
(|IT|) and p3|) are called Grassmann variables. Thus el- 
ements of the Grassmann algebra can be second order 
polynomials at the most. 



7(77,77*) = A + 57/ + C77* + 157777* 



(15) 



and the complex conjugate of a product of two elements 
is equal to the product of the conjugates written in the 
reverse order. 
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Using Eqs. (|TTj) and (|I3|) . we can write the overlap 
between the two coherent states as 



This imphes that 



(77h) = 1 +7^*77 = 6"' (16) 

Integration of the Grassmann variables is defined by, 

(17) 
(18) 

The differential elements dt] and drj* anticommute with 
each other. Using Eqs. jlll , Hll , jlll and JTH]) it is 
straightforward to show that 



drj — J drf — 
drirj = / drfrf = 1. 



dr]d'n*e^''^\7^){'q\ = 1 



(19) 



which is the closure relation for coherent states. 

Differentials of the Grassmann variables are defined as, 

^/(77, V*)^B + D7f, -^firj, rj*)^C- i5r/.(I10) 



d d 
dr] dr)* 



d d 
drj* dr] 



(111) 



Taking integral of f{ri,ri*) with respect to rj or rj* and 
comparing with Eqs. pil|) . we obtain the operator iden- 
tities 



(112) 



Using Eqs. dlZl), jlSj) and (iI12|) . it is straightforward to 
see that for any N x N matrix A, 



J I?(7/*77)e^'^'''^'^''^ =Det[^] 



(113) 



where T>{r]*ri) — JJ^ dT]*d7]i 
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